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SECTION  1 


INTRODUCTION 

In  this  report,  we  develop  a  solution  to  the  axisymmetric  problem  of  a 
tunnel  in  a  Mohr-Coulomb  medium  with  internal  and  far-field  applied  pres¬ 
sures.  The  solution  presented  herein  extends  existing  approaches  by  address¬ 
ing  generalized  stress  conditions  and  incorporating  an  arbitrary  dilatancy 
parameter.  This  parameter  may  be  varied  between  certain  limits,  but  is  as¬ 
sumed  to  be  a  constant  for  any  given  problem. 

1.1  BACKGROUND. 

The  basic  approach  to  this  problem  was  defined  by  Newmark  (Reference  1) 
prior  to  the  widespread  availability  of  high-speed  digital  computers  and  in 
response  to  the  need  for  a  method  for  designing  underground  structures  to 
resist  the  effects  of  nuclear  weapons.  In  the  development  of  his  solution. 
Newmark  assumed  an  elastic-plastic  medium  with  volume  constancy  in  a  plane 
strain  configuration.  Newmark 's  work  was  expanded  by  Hendron  and  Aiyer 
(Reference  2),  who  provided  the  solution  for  a  dilatant  material  which  obeys 
an  associated  flow  rule.  Newmark,  as  did  Hendron  and  Aiyer.  assumed  that  the 
out-of-plane  stress  (Oz)  was  the  intermediate  principal  stress  and  that  the 
circumferential  stress  was  the  minimum  (most  compressive)  principal  stress. 
The  out-of-plane  stress  assumption  is  valid  for  a  limited  range  of  material 
properties  and  loads;  however,  it  is  reasonable  for  many  tunnel  situations 
and  yields  a  relatively  simple  solution  for  the  stress  and  displacement 
fields  because  only  one  plastic  zone  develops. 

Florence  and  Schwer  (References  3  and  4)  generalized  the  Hendron  and 
Aiyer  solution  by  eliminating  the  requirement  that  <jz  be  the  intermediate 
principal  stress.  Their  objective  was  to  provide  an  analytical  solution  for 
use  in  validating  finite  element  calculations  in  support  of  experimental  and 
theoretical  investigations  of  buried  cylinders. 

Eliminating  the  requirement  that  Oz  be  the  intermediate  principal 
stress  substantially  increases  the  complexity  of  the  solution  by  requiring 
one  to  deal  with  multiple  yield  conditions  in  a  single  problem.  As  stated  in 
Reference  3,  these  "depend  on  the  relative  magnitudes  of  the  principal 
stresses,  and  these  relative  magnitudes  depend  on  the  values  of  Poisson's 
ratio  and  the  friction  angle."  Florence  and  Schwer  developed  the  solutions 
for  two  different  cases.  For  Case  I  (Reference  3),  there  are  three 
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different  plastic  zones  (each  with  a  different  set  of  governing  equations)  in 
the  material  surrounding  the  tunnel.  The  problem  geometry  is  illustrated  in 
Figure  1  (from  Reference  3).  For  Case  II  (Reference  4),  there  are  two 
different  plastic  zones  (each  with  a  different  set  of  governing  equations)  in 
the  material  surrounding  the  tunnel.  Case  II  geometry  is  shown  in  Figure  2. 

Detournay,  St.  John,  and  Van  Dillen  (Reference  5)  generalized  the 
Hendron  and  Aiyer  solution  to  allow  the  use  of  an  arbitrary  dilatancy  parame¬ 
ter  but  their  solution  maintained  the  requirement  that  Oz  be  the  intermediate 
principal  stress.  A  similar  solution  was  developed  independently  by  Merkle 
(Reference  6),  although  we  have  not  found  it  in  published  form.  Use  of  the 
arbitrary  dilatancy  parameter  allows  the  material  to  undergo  a  volume  change 
that  is  intermediate  between  the  volume  constancy  imposed  by  the  Newmark 
solution  and  the  strictly  associated  flow  rule  dilatancy  imposed  by  Hendron 
and  Aiyer.  Detournay  and  St.  John  (Reference  7)  assert  that  the  dilatancy 
parameter  should  be  a  variable  that  is  a  function  of  accumulated  plastic 
shear  strain  (and  possibly  of  mean  pressure).  However,  the  use  of  such  a 
variable  would  require  the  application  of  numerical  integration  in  obtaining 
a  solution. 

1.2  SCOPE  AND  ORGANIZATION  OF  THE  REPORT. 

The  solutions  discussed  in  the  previous  section  are  “exact"  only  for 
axi symmetrical ly  applied  static  loads.  Nevertheless,  these  methods 
(particularly  those  of  Newmark  and  of  Hendron  and  Aiyer)  have  proven  very 
useful  in  investigating  the  effects  of  nuclear  weapons  on  underground 
structures  constructed  in  rock.  They  have  been  used  extensively  by  members 
of  the  DNA  community  for  developing  preliminary  designs  of  tunnel  support 
systems,  preliminary  hardness  assessments  of  existing  structures,  and 
estimates  of  targeting  requirements. 

Because  they  may  be  programmed  on  a  personal  computer  and  then  easily 
and  rapidly  applied,  we  expect  these  "exact"  solutions  will  continue  to  be 
widely  used  in  the  future  to  obtain  preliminary  structural  response  estimates 
prior  to  performing  detailed  finite  element  calculations.  Consequently,  we 
have  extended  the  "exact"  solution  one  step  further  by  combining  the  general¬ 
ized  stress  state  of  References  3  and  4  with  the  concept  of  arbitrary  dila¬ 
tancy  from  References  5,  6,  and  7. 

The  solution  is  presented  in  the  format  used  by  Florence  and  Schwer  in 
References  3  and  4,  insofar  as  possible.  We  have  used  the  same  notation,  and 
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equations  are  presented  in  the  same  formats  (and,  to  the  extent  possible, 
with  the  same  numbers)  as  used  by  Florence  and  Schwer.  Many  of  our  equations 
are  identical  to  theirs  and  much  of  the  introductory  material  is  largely  a 
review  of  material  in  Reference  3.  Each  symbol  is  defined  where  first  used, 
and  definitions  of  the  more  frequently  used  symbols  are  summarized  in  Table 
1.  We  have  adopted  the  Florence  and  Schwer  sign  convention  by  defining 
tensile  stresses  and  strains,  applied  pressures,  and  unconfined  compressive 
strengths  as  positive.  This  approach  allows  the  reader  to  compare  equations 
with  and  without  arbitrary  dilatancy  throughout  the  derivation. 

In  general,  our  approach  is  to  replace  the  friction  parameter  *N" 
(defined  in  Equation  2  below)  by  the  arbitrary  dilatancy  parameter  "M"  in 
expressions  involving  plastic  strains  or  displacements.  It  can  be  demon¬ 
strated  that  for  associated  flow  (M  =  N),  the  present  solution  is  identical 
to  that  of  References  3  and  4. 

The  general  solution  for  initial  yield  and  the  conditional  inequalities 
leading  to  Cases  I  and  II  are  presented  in  Section  2.  The  solutions  for 
Cases  I  and  II  are  developed  in  Sections  3  and  4,  respectively.  Step-by-step 
procedures  for  applying  the  closed-form  solutions  and  example  problems  that 
demonstrate  their  use  are  provided  in  Section  5.  With  the  permission  of  the 
publisher.  Reference  3  is  presented  in  its  entirety  in  Appendix  A.  It  is  our 
understanding  that  the  Case  II  solution  of  Reference  4  was  developed  under  a 
DNA  contract,  but  was  not  published.  That  reference  is  provided  in  Appendix 
B  with  the  concurrence  of  the  senior  author. 
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Table  1.  Definitions  of  frequently  used  symbols. 

a  -  interior  radius  of  the  opening 
E  -  Young's  modulus  for  medium 

G  -  shear  modulus  for  medium 

M  -  arbitrary  dilatancy  parameter 

N  -  friction  parameter 

Pa  -  internal  pressure  acting  on  hole  boundary 

p*  -  negative  of  the  radial  stress  at  the  elastic-plastic  boundary 

Pb  -  far-field  pressure  (compressive  stress)  at  large  radius 
Pb  _  far-field  pressure  at  initial  yield 

Pb  -  far-field  pressure  at  which  R  =  R 

Pb  -  far-field  pressure  when  the  inner  plastic  zone  begins  to  form 

under  Case  II  yield  conditions 
r  -  arbitrary  radius  to  any  point  in  the  medium 

R  -  radius  to  the  elastic-plastic  boundary 

R  -  maximum  radius  at  which  the  radial  and  out-of-plane  stresses 
are  equal 

R  -  minimum  radius  at  which  the  radial  and  out-of-plane  stresses 
are  equal 

R'  -  radius  to  the  elastic-plastic  boundary  when  Pb=P^ 

Ou  -  unconfined  compressive  strength  of  medium 

Or  -  stress  in  radial  direction 

O0  -  stress  in  circumferential  direction 

Oz  -  stress  in  out-of-plane  direction 

Er  -  strain  in  radial  direction 

£0  -  strain  in  circumferential  direction 

Ez  -  strain  in  out-of -plane  direction 

Note:  Strains  are  further  distinguished  by  the  superscripts  (e) 

to  denote  elastic  component  and  (p)  to  denote  plastic  component 
V  -  Poisson's  ratio  for  medium  (0  <  v  <  0.5) 

(p  -  friction  angle  for  medium 
X  -  flow  constant  of  proportionality 
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SECTION  2 


GENERAL  CONDITIONS 

In  this  section,  the  general  solution  is  developed  for  initial 
yielding  of  the  medium  around  the  tunnel.  The  conditional  inequalities  that 
lead  to  Cases  I  and  II  are  identified.  All  equations  in  this  section  are 
identical  to  those  in  the  corresponding  portion  of  Reference  3. 

2.1  INITIAL  LOADING-OUTER  ELASTIC  ZONE. 

Initially,  let  the  internal  and  far-field  (theoretically  at  infinity) 
pressures  pa  and  Pb  be  equal.  If  these  pressures  are  increased  simultane¬ 
ously  until  the  material  yields,  initial  yielding  will  be  governed  by  the 
yield  function 


f  =  Oe  -  NOz  +  0^=0 


(1) 


where 

N=LLiin^ 

l-sin(p 


(2) 


Substitution  of  stresses  (Or  =  ae  =  -Pb.  Oz  =  -2vpb)  in  Equation  1  gives 


Pa=Pb  = 


l-ZNv 


(3) 


And  if  Equation  3  is  satisfied,  the  material  will  yield  completely  if  2Nv  <  1 
(such  yielding  is  precluded  if  2Nv  >1);  consequently,  to  insure  elastic 
conditions  in  the  far-field,  the  internal  pressure  will  be  limited  such  that 

Pa  <  — —  when  2Nv<l  (4) 

"  1-2NV 


though  no  such  limitation  is  required  if  2Nv  >  1. 

For  the  more  general  condition,  once  loading  with  pa  =  Pb  reaches  the 
desired  value  of  Pa  (e.g.,  the  yield  strength  of  a  tunnel  liner  or  the  crush 
strength  of  backpacking  material),  pb  can  be  increased  to  some  higher  value, 
Pb  >  Pa-  In  the  elastic  condition,  stresses  and  strains  are  given  by 

a2  a2 

=-Pb-^{Pb-Pa)-7’  Oe=-Pb-(Pb-Pa)  — -  O,  = -2vPb  (5) 
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e,=0  (6) 


ZGe^  =-(l-2v)Pb  +  (Pb-Pa)\,  2Gee  =  -(l-2v)Pb-(Pb-Pa)-^. 

ft  pt 

Note  that  the  negative  sign  in  the  third  of  Equations  5  was  omitted  in 
Reference  3. 

2.2  INITIAL  YIELDING. 

Yielding  will  always  begin  at  the  inner  surface,  that  is,  at  r  =  a. 
The  far-field  pressure  associated  with  initial  yield  is  defined  as  Pj,.  In 
Equations  5,  Or  >  -Pb.  oe  <  -pb.  and  <s^  >  -pb.  Thus,  the  minimum  (most 
compressive)  principal  stress  will  always  be  in  the  circumferential  direc¬ 
tion,  but  the  maximum  principal  stress  may  be  either  radial  or  out-of-plane. 
If  the  internal  pressure  (and  hence  the  radial  stress  at  the  opening)  is 
lower  in  magnitude  than  2vPb  (the  out-of-plane  stress),  then  ae  <  Oz  <  Or  at 
the  onset  of  yield  (i.e.,  Or  is  the  maximum  principal  stress),  so  the  corres¬ 
ponding  yield  function  is 

f  =  00  -NOr  +au  =  0  (7) 

The  radial  and  circumferential  stresses  at  the  interior  are  computed  at  yield 
from  Equations  5.  Substitution  of  those  expressions  in  Equation  7  gives  a 
value  for  the  external  pressure  P^  at  the  onset  of  yield. 

Pb  =^[(N  +  l)Pa-^Ou]  (8) 

In  order  to  satisfy  the  inequality  pa  <  2vPb,  it  is  required  that 

Pa< — ; — (N+1)V<1  (9) 

®  l-(N  +  l)v  '  ' 


Alternatively,  where  Equation  1  governs  the  yield  condition,  the  yield 
pressure  is  given  by 


Pb  = 


Pa^^u 
2(1 -Nv) 


Nv<  1 


(10) 


For  Nv  >  1.  the  yield  condition  of  Equation  7  rather  than  that  of  Equation  1 
applies;  consequently.  Equation  8  is  applicable.  With  Equation  10,  the  fol- 


& 


Governing  conditions  for  the  above  situations  are  summarized  in  Table  2. 
which  is  equivalent  to  Table  1  in  both  References  3  and  4. 
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Table  2.  Conditions  at  initial  yield. 


Yielding  Case  Stress  Order  Property 
at  r  =  a  Relations 


Internal 

Pressure 


Far-field 
Pressure  at 
Initial  Yield 


Case  I 


P.  <  2vp, 


Co<o  <a,  (N-^l)v<l 


l-(N  +  l)v 


1<  Nv 

Nv<l<{N  +  l)v  p.>0 


Case  II 


2vp,<p.<p,  0e<0,  <0,  j<Nv.  (N  +  l)v<l 


_  p  - 

(N+l)v  "  2(1 -Nv) 


0<  Nv<  I 


- y — <  p  <  — y — 

l-(N  +  l)v  *  1-2NV 
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SECTION  3 


CASE  I  SOLUTION 

As  stated  earlier,  in  Case  I,  three  separate  plastic  zones  (with 
differing  sets  of  governing  equations)  form  in  the  medium  surrounding  the 
tunnel.  The  governing  equations  are  developed  in  the  four  subsections  that 
follow.  In  Section  3.1,  the  equations  associated  with  initial  growth  of  the 
plastic  zone  (resulting  from  increasing  pb)  are  developed.  This  is  followed 
by  three  subsections  in  which  the  conditions  governing  the  formation  of  the 
outer,  intermediate,  and  inner  plastic  zones  are  developed. 

3.1  INITIAL  GROWTH  OF  THE  PLASTIC  ZONE. 

For  Case  I,  the  yield  condition  of  Equation  7  controls  in  the  outer 
plastic  zone.  From  equilibrium  considerations  we  write 

r-^  +  CTr-a0  =  O  (13) 

dr 


Circumferential  stress  can  be  determined  from  Equation  7  and  substituted  in 
Equation  13.  The  result  is  then  solved  such  that  the  radial  stress  at  the 
tunnel  surface  (r  =  a)  is  equal  to  -pa 


N-1 


a<  r  <  R 


(14) 


The  flow  rule  for  the  yield  condition  of  Equation  7  is  given  by 


p(p)  - 1 


®  aa, 


e 


(15) 


Note  that  the  first  relationship  in  Equations  15  constitutes  the  first  depar¬ 
ture  from  the  derivation  of  Reference  3.  It  differs  from  the  corresponding 
equation  in  Reference  3  only  in  the  substitution  of  M  for  N  in  the  radial 
strain  component.  The  arbitrary  dilatancy  factor  M,  introduced  in  References 
5  and  6,  can  assume  any  value  between  one  and  N.  Dots  are  used  here  to 
designate  differentiation  with  respect  to  external  pressure.  From  this  flow 
ru^e,  the  following  strain  equations  may  be  written  (gi  and  g2  are  unspecified 
functions  of  r). 


=  0, 


e[.P’>MEj,‘’^  =  gi(r).  e[P>=g2(r)  (16) 
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At  the  elastic-plastic  boundary  (r  =  R),  there  can  be  no  plastic  strains,  so 
we  can  write 


e<P^Me^P^  =  0,  e^P^=0  (17) 


Since  strain  increments  are  the  sums  of  their  elastic  and  plastic  components, 
we  can  write,  for  strain  increments  and  total  strains,  respectively. 


p  =  p(®)  X  pip) 

Ee  Eq  +Eq  , 

£2  -  e^®^  +  4’’^ 

(18) 

Er  =  4®^  + 

c  =  pi®)  ♦  piP) 

*^0  E0  ♦Cq  , 

E^  =  E^®^  +  E^'*^ 

(19) 

Now. 

from 

i  .e. 

plane  strain 
Equations  17 

requires  that  the  total 
we  can  conclude  that  the 

out-of-plane  strain  be  zero,  so 
elastic  component  is  also  zero. 

e‘P)=0 

(20) 

Elastic  strains  are  related  to  stresses  by  Hooke's  law. 


EE<,®>  =  Op  -v(ae+oJ 

(21) 

Ee^®’  =a0-v(o^  +  ap) 

(22) 

Ee^®>  =a^-v(ar  +09) 

(23) 

And.  since  Cz  is  zero  in  the  elastic  region  (outside  the  elastic-plastic 
boundary  at  r  =  R).  from  Equation  23 


o,  =  v(a,  +O3) 


(24) 


From  the  yield  condition  of  Equation  7.  circumferential  stress  at  the 
elastic-plastic  boundary  can  be  expressed  as  a  function  of  the  radial  stress. 
Using  this  and  Equation  24,  Equations  21  and  22  can  be  used  to  express  the 
radial  and  circumferential  elastic  strains  as  functions  of  radial  stress: 

2G4®>=[l-(N  +  l)v]a, +VO,  (25) 
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2Ge^®^  =[N-(N  +  l)v]ar  -(l-v)au 


(26) 


Note  that  for  conputational  convenience,  the  shear  modulus  is  substituted  for 
Young's  modulus  at  this  point,  as  allowed  by  the  relationship  E=2G(l+v),  and 
strain  terms  in  subsequent  equations  are  frequently  expressed  in  the  form  of 
the  left-hand  side  of  Equation  26. 

For  the  axisymmetric  configuration,  the  following  strain-displacement 
relationships  apply 


Er 


(27) 


from  which  the  compatibility  equation  can  be  written  as 


r 


dr 


+  £0-6^ 


=  0 


(28) 


Using  Equations  19  to  separate  strains  into  elastic  and  plastic  compo¬ 
nents,  Equations  17  to  express  radial  strain  in  terms  of  circumferential 
strain,  and  Equations  25  and  26  to  substitute  for  the  elastic  strains,  we  can 
rewrite  Equation  28  as 


de 


(P) 


dr 


(M  +  l)e 


(P) 

e 


1-v 

2G 


(Nn)[(N-l)p3  +  oJ 


V  a 


N-l 


(29) 


By  inspection,  it  can  be  seen  that  this  differs  from  the  corresponding 
expression  in  Reference  3  only  in  the  substitution  of  (M+1)  for  (N+1)  on  the 
left-hand  side,  and  indeed  degenerates  to  the  identical  equation  for  the  case 
where  M  =  N.  The  solution  to  Equation  29  is  readily  determined  by  direct 
integration  over  the  interval  from  r  =  r  to  r  =  R.  The  resultant  solution  is 


1-v 

2G 


N>1 
M  +  N 


[(N-l)Pj+oJ 


(-] 


N-l 


(30) 


The  counterpart  of  equation  30  in  Reference  3  contains  a  typographical  error 

N  +  1 

in  that  the  second  term  on  the  right-hand  side  should  read  .  The  plastic 

2N 

component  of  radial  strain  can  be  obtained  from  Equations  17  although  at  this 
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point,  R  remains  unknown. 


If  the  material  inside  of  radius  r  =  R  is  assumed  to  be  replaced  with 
a  pressure  p,  it  follows  from  Equations  5  that 


Or(R)  =  -P. 


Oe(R)  =  -2Pb  +  P. 


o^{R)  =  -2vpb 


(31) 


Since  the  yield  condition  of  Equation  7  is  applicable  at  r  =  R, 


Pb  =  2t(N  +  l)p  +  Ou] 


gr(R)  =  --' 

N  +  1 


(32) 

(33) 


Equating  the  radial  stress  from  Equation  33  to  the  radial  stress  computed 
from  Equation  14  for  r  =  R  yields 


f-1 

u> 


N-1 


_  2  (N-l)Pb+o, 


(N-l)Pa+a, 


(34) 


When  (N+l)v  >  1,  Oz  will  be  the  intermediate  principal  stress  through¬ 
out  the  plastic  zone  and  the  foregoing  procedure  is  adequate  to  completely 
define  the  solution.  This  corresponds  to  the  Hendron  and  Aiyer  solution  as 
modified  by  Detournay  and  St.  John  to  include  arbitrary  dilatancy. 

If  (N+l)v  <  1,  the  elastic-plastic  boundary  still  forms  at  r  =  R_  but 
this  condition  exists  only  as  long  as  Oz  is  the  intermediate  principal 
stress,  i.e.  as  long  as  Or  <  -vau/[l-(N+l )v].  From  Equation  14.  we  see  that 
Or  is  maximum  (least  compressive)  at  r  =  a  and  minimum  (most  compressive)  at 
r  =  R.  This  implies  that  as  the  loading  increases,  a  plastic  zone  forms,  in 
which  the  radial  and  out-of-plane  stresses  are  equal.  In  Figure  1,  this  is 
the  annular  region  between  r  =  R  and  r  =  R,  where  R  and  R  are  defined  as  the 
minimum  and  maximum,  respectively,  radii  where  Or  =  Oz.  At  r  =  R  then,  we 
have 


OrfR)  =  a^(R)  = - ^ 

’  1-(N  + 


(N+l)v 


(35) 


and  from  plane  strain  considerations,  the  circumferential  stress  is 

aJR)  =  -JLl^ 

'  l-(N+l)v 


(36) 
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We  compute  the  external  pressure  at  which  R  =  R  from  Equations  31  and  35  as 


Pb= - - 

2[l-(N  +  l)v] 

which  can  then  be  substituted  in  Equation  34  to  compute  R 
/  -  \N-i 

(i-2v)qu 

,aj  [l>(N  +  l)v][(N-l)p,+aJ 


(37) 


(38) 


Note  that  Equation  38  does  not  involve  Pb^  from  which  we  can  conclude 
that  as  pb  is  increased  beyond  Pb,  R  remains  constant. 

3.2  INCREASED  LOADING-OUTER  PLASTIC  ZONE. 

As  discussed  earlier.  aa<o..<a,  at  r  =  R.  Therefore  there  must  be  a 
third  plastic  zone  between  r  =  R  and  r  =  R,  as  shown  in  Figure  1,  which 
develops  when  Pb  is  increased  beyond  Pb.  In  this  region,  the  yield  condition 
of  Equation  1  applies  and  the  flow  rule  (with  arbitrary  dilatancy)  is 

4'’^  =  0,  4'*^  = -MX  (39) 

In  the  same  manner  that  strain  relations  were  developed  from  the  flow  rule  of 
Equations  15,  strain  relations  can  be  derived  from  Equations  39  as 

+  Me^^^  =  0,  =  0  (40) 

and  from  the  plane-strain  condition, 

6^=4"^-  ee=4®^  +  ^4"’’  E,  =0  (41) 

Also,  in  the  manner  that  Equations  25  and  26  were  developed  using  the  Hooke's 
law  relationships  (Equations  21  through  23).  Equations  41  can  be  used  with 
Equation  1  and  Hooke's  law  to  obtain  expressions  for  radial  and  circumferen¬ 
tial  stress  in  the  outer  plastic  zone 

[MN  +  l-(M  +  l)(N+l)v]o,  =[MN  +  l-(M  +  N)v]2Ge^  +M(N  +  l)v2Gee+(M-l)va,  (42) 
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[MN  + 1  -  (M  +  1)(N  + 1) v]ae  =  N(M  +  l)v2Ge,  +  MNZGEg  -  [1  -  (M  +  l)v]a. 


(43) 


and  with  Equations  13  (equilibrium)  and  27  (strain-displacement),  this 
results  in  the  following  differential  equation  for  radial  displacement 


^  ,  r  5!i[l .  -P'u  = 

I  ?  I  kill  r  J  r  .j4| 


2Gu  =  Ajr^‘  +  A^r  +  A-  r  (46) 

The  corresponding  strains  are 

2G—  =  2Ge,  =  " YzA^r'^*''  +  A  (47) 

dr 

2G-  =  2Ge  =A,r^‘''+/Lr'^'''+A  (48) 

r 


where 


and  Ai  and  ^2  3re  constants,  yet  to  be  determined. 

Compared  to  the  associated-flow-rule  formulations  of  Reference  3. 
Equations  46,  47,  and  48  involve  additional  factors  and  non-identical  powers 
of  r.  In  spite  of  these  complexities,  for  the  case  where  M  =  N,  Yi  =  Y2  =  (5, 
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and  p  is  identical  to  the  term  a  used  by  Florence  and  Schwer.  Consequently, 
Equations  46,  47,  and  48  become  identical  to  their  Reference  3  counterparts. 

We  can  employ  the  equilibrium  and  compatibility  conditions  at  the 
elastic-plastic  boundary  (r  =  R)  to  equate  radial  stresses  and  circumferen¬ 
tial  strains  across  the  boundary.  From  Equations  1  and  5,  the  elastic-zone 
radial  stress  is 

o,{R)  =  -2(l-Nv)Pi,+a,  (50) 

from  which  the  circumferential  strain  can  be  obtained,  using  Equation  26,  as 

2Gee(R)  =  -2(N-l)vPb-Ou  (51) 

Invoking  compatibility  of  radial  displacement  (or  circumferential  strain),  we 
can  equate  the  right-hand  sides  of  Equations  48  and  51,  resulting  in 

[o.-a-2Nv)p,]  (52) 

From  equilibrium,  the  elastic  radial  stress  from  Equation  50  is  equated  to 
the  plastic  zone  stress  from  Equation  42,  and  by  substituting  strains  from 
Equations  47  and  48 

bjAjR^-'^  +  b^A^R'^'’' =0^  (53) 


where 

bi  =  M(N  + 1)  V  +  [MN  + 1  -  (M  +  Nlvjy^, 
b2  =  M(N  +  l)v-[MN+l-(M  +  N)vl72. 

Cj  =  2^^^[MN+l-(M  +  l)(N+l)v][o„  -(l-2Nv)pJ 


Since  Equations  52  and  53  are  linear  in  Ai  and  Az.  we  can  solve  simultaneously 
(in  terms  of  unknown  R)  as 
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At  the  assumed  plastic-plastic  boundary  at  r  =  R  and  in  the  plastic 
zone  between  r  =  R  and  r  =  R,  the  radial  and  out-of-plane  stresses  are  equal 
and  the  yield  conditions  of  Equations  1  and  7  are  both  valid.  Using  Equation 
7,  we  can  now  write 

Nar(R)-Oe(R)  =  au  (54) 

Since  the  stresses  from  Equations  42  and  43  are  valid  at  r  =  R.  we  can 
substitute  them  in  Equation  54  to  get 

b,A,R^i‘Wb,A,R‘^2-i  =  c  (55) 

3  1  2  n-2Nv  “ 


where 


b3  =  M[(N  +  l)v  -1]  +  [MN  + 1  -  "’ri  -  N  -  1)v]Yi 
b4  =  M[(N  +  i  jv  -1]  -  [MN  + 1  -  (2M  +  N  +  1)v]y2 
C3  =[MN+l-(M+l)(N+l)v] 


The  expressions  for  Ai  and  A2  in  Equation  53.1  can  be  substituted  in  Equation 
55  to  define  the  ratio  of  R  to  R. 


bsBi 


l-2v 

n-2Nv  “ 


(56) 


Since  the  ratio  cannot  be  solved  for  explicitly,  an  implicit  scheme  will  be 
required,  but  rapid  convergence  is  expected.  The  radius  R  remains  unknown 
npnding  consideration  of  the  region  from  r  =  R  to  r  =  R. 

3.3  INCREASED  LOADING— MIDDLE  PLASTIC  ZONE. 

As  noted  above,  both  of  the  yield  conditions  (Equations  1  and  7)  are 
valid  in  the  region  R  <  r  <  R,  and  are  restated  here  to  continue  the  paral¬ 
lelism  with  Reference  3. 

f;  =  09-No, +0,  =0  (57) 
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fj  =ae  -No,  +a„  =0 


(58) 


From  Equation  57  we  obtain  an  expression  for  circumferential  stress  to 
substitute  into  the  equilibrium  equation  (Equation  13).  The  resulting 
differential  equation  is 

r^  +  a,(l-N)  =  -o„  (58.1) 

dr 

for  which  the  solution  is 


(59) 


If  Or (R)  and  R(from  Equations  35  and  38)  are  substituted  in  Equation  59,  it 
becomes 


Or 


N-1 


Va 


N-l 


(60) 


Since  Equations  57  through  60  are  independent  of  the  flow  rule,  they  are 
identical  to  their  counterparts  in  Reference  3. 


We  recall  that  Equation  42  (radial  stress)  is  valid  in  the  region  R  <r 
<  R  and  Equation  60  is  valid  for  R  <  r  <  R.  Since  we  must  assure  continuity 
of  radial  stress  at  r  =  R,  we  can  equate  the  two  expressions  (evaluated  at  r 
=  R)  to  obtain  an  expression  for  R 


biAiR^>'^  +  b2A2R‘^2'^  =C3 


N(l-2v) 

(N-l)(l-2Nv) 


Ou  "^3 


(61) 


with  bi.  bi,  and  C3  as  defined  in  Equations  53  and  55.  The  expressions  for  Ai 
and  A2  from  Equations  53.1  can  be  substituted  in  Equation  61: 


biBi 


+  5^62 


=  Cl 


N(l-2v) 

(N-l)(l-2Nv) 


Ou  -C; 


Pa  + 


N-l. 


(62) 


Since  R/R  was  determined  in  Equation  56.  R  can  be  computed  from  Equation 
62.  R  can  then  be  obtained  from  R/R.  Then,  Aj  and  A2  can  be  calculated 
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explicitly  from  Equations  53.1,  which  is  repeated  here  as  Equation  63.  again 
to  continue  the  parallelism  with  Reference  3. 


c^-b^qA  1  ^__B, 


^  b,-b, 


R-Tf.-I 


_  B3 


>-Y2-1 


(63) 


Strains  and  displacements  in  the  region  R  <  r<  R  can  then  be  computed 
from  Equations  46  through  48;  stresses  are  available  from  Equations  1.  42, 
and  43. 


For  the  yield  conditions  of  Equations  57  and  58,  the  flow  rule  with  ar¬ 
bitrary  dilatancy  is 


pip)  _  y  .  y  ^2  -  uY 


(64) 


p(p)  _  j  ^ 


(65) 


(66) 


from  which 


g(r )  =  e|f  *  +  Meg*^  +  =  0 


(67) 


We  designate  the  plastic  strains  as  ef\  Eg’^,  when  r  =  R,  so  that 
g(r)  =  +  MCg’*  + 


(68) 


Since  Equation  40  is  applicable  in  the  region  R<  r  <  R,  g(r)  =  0,  from  which 
gfP)  +  Me^gP^  +  =  0  (69) 

and  with  e[®^  = from  the  plane  strain  condition. 


Ep  =  -Mcg  +  Me^®*  +  Ep®’  +  E^®* 


(70) 
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With  Or  =  Oz  given  by  Equation  60  and  using  the  yield  condition  of 
Equation  57  we  can  write  Hooke's  law  from  Equations  21  through  23  as 


Ee'*>=Ee^*^  =  [l-(N+l)v]q, +  vq. 

(71) 

Ee^®^  =(N-2v)qr 

(72) 

Elastic  strains  from  Equations  71  and  72  can  be  substituted  into  Equation  70. 
Substitution  of  the  radial  strain  expression  thus  obtained  into  the  compati¬ 
bility  equation  (Equation  28)  yields  a  differential  equation  in  terms  of  the 
circumferential  strain 


_d_ 
dr  '■ 


2G— fr-e 


(N-l)(l  +  v) 


M*N-1 


^(M>2)a-2v)o,^H 

(N-l)(l  +  v) 


(73) 


This  differential  equation  is  directly  integrable,  and  integration  over  the 
limits  of  r  to  R  results  in  the  following  expression  for  ee 


2Gee  =  2Gee(R)- 


+ 


MN-^2-2(M-^N4l)v 
(M  +  N)(N-1)(1  +  V) 


[(N-l)p.+aJ 


R 

r 


(M^2)(l-2v)q„ 
(N-l)(M  +  l)(l  +  v) 


(74) 


e0(  R )  can  be  obtained  from  Equation  48  since  we  now  have  values  for  Ai  and 
A2.  The  radial  strain  is  obtained  from  Equations  70  through  72  giving 


2Ge^  = -2GMeg  + 


MN  +  2-2(M  +  N+l)v  M-2v 


1  +  V 


q,  — 


1  +  V 


(75) 


3.4  INCREASED  LOADING— INNER  PLASTIC  ZONE. 

For  an  applied  pressure  in  excess  of  Pb,  stresses  in  the  inner  plastic 
zone,  a  <  r  <  R,  remain  constant.  Radial  stress  is  given  by  Equation  14  and 
the  radial  and  out-of-plane  plastic  strain  cotrponents  are  given  by  Equations 
17,  Equation  29  is  valid  in  this  region,  but  we  must  perform  the  integration 
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from  r  to  R  f or  this  condition,  which  yields  the  circumferential  plastic 
strain  as 


2Ge§" 


^N-l 


/  - 


(76) 


Elastic  strains  are  available  from  Equations  20  through  26.  In  order  to 
compute  the  circumferential  strain  at  r  =  R  f rom  Equation  76,  we  can  compute 
the  total  strain  from  Equation  74  evaluated  at  r  =  R  and  subtract  the  elastic 
component.  Circumferential  and  out-of-plane  stresses  are  determined  using 
Equations  7  and  24. 
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SECTION  4 


CASE  II  SOLUTION 

As  stated  earlier.  Florence  and  Schwer  presented  the  solution  for  Case 
II  in  Reference  4.  Since  that  reference  is  a  "stand-alone"  paper,  they  re¬ 
peated  some,  but  not  all.  of  the  equations  from  Reference  3  with  different 
equation  numbers.  Consequently,  from  this  point  forward,  we  do  not  use  the 
Florence  and  Schwer  equation  numbers,  but  we  continue  to  follow  their  format. 
For  the  convenience  of  the  reader,  several  of  the  equations  developed  earlier 
are  repeated  in  this  section. 

As  discussed  in  both  References  3  and  4,  there  are  at  most  two  plastic 
zones  for  Case  II.  The  problem  geometry  is  shown  in  Figure  2  (from  Reference 
4). 


4.1  INITIAL  GROVfTH  OF  THE  PLASTIC  ZONE. 

As  shown  in  Table  2.  the  order  of  the  principal  stresses  at  initial 
yield  is  ae  <  Or  <  Oz  for  Case  II.  Therefore,  the  yield  condition  of 
Equation  1  applies. 


f  =  aQ-Naz  +  au  =0  (1) 

For  this  condition,  with  arbitrary  dilatancy,  the  flow  rule  was  given  in 
Equations  39, 

4^>  =  0.  4”^  = -MX  (39) 

from  which  the  plastic  strain  relationships,  given  by  Equations  40, 

gjp)  +  =  0,  =  0  (40) 

were  derived.  Incremental  and  total  strains  are  given  by  Equations  18  and 
19. 


Er  =£r^£r 


(p) 


P  _p(e)  p(p) 


E,  =  E^^^  +  E^'^^ 


(18) 


£  =  £(®)  +  £(P) 

tp  Cp  tp  , 


p  _p(e)  (p) 

^0  *^0  *^0 


E^  =  E^®*  +  E^'’' 


(19) 
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from  which  the  three  components  of  total  strain,  given  by  Equations  41, 


er=4®^.  (^1) 

M 

are  obtained.  The  yield  condition  of  Equation  1  allowed  us  to  eliminate  Oz 
from  the  Hooke's  law  expressions  of  Equations  21,  22,  and  23, 


Eej.®^  =  ar  -\{aQ*aj) 

(21) 

E4®^  =  a0-v(a^  +  ap) 

(22) 

Eei®^  =  a^-v(Op+a0) 

(23) 

which  led  to  the  expressions  for  radial  and  circumferential  stress  in  terms 
of  strain  in  Equations  42  and  43. 


[MN  + 1  -  (M  +  1)(N  +  l)v]ar  =  [MN  + 1  -  (M  +  N)v]2Ger  +  M(N  +  l)v2Gee  +  (M  -l)vOu 


(42) 


[MN  + 1  -  (M  +  1)(N  + 1) v]a0  =  N(M  +  l)v2Ger  +  MN2Ge0  -  [1  -  (M  +  l)v]au 


(43) 


Equation  13. 

dOp 

r — - 
dr 


+  -O0  =0 


(13) 


(equilibrium)  and  Equations  27, 
du 


Er  = 


dr 


u 


(27) 


(strain-displacement)  were  then  used  to  obtain  the  differential  equation  of 
displacement  given  by 


2  d^u  ,du,,  .  (H-N)v„2,  2,._  (l-2v)Ou 


+  r— [1 
dr'^  dr 


MN 


P1-P^^=~ 


MN  2G 


(44) 


The  solution  of  Equation  44  is 
2Gu  =  Ajr^'  +  +  A-  r 
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(46) 


and  the  associated  strain  relationships  are 


2G—  =  2Ge,  =  +  A 

dr 

(47) 

2G^  =  ZGEe  =  A,r^'‘‘  +  A^r'^*''  +  A 

(48) 

Continuity  of  displacement  (or  circumferential  strain)  and  radial  stress  at 
the  elastic-plastic  boundary  (r  =  R)  then  allowed  us  to  write  Equations  52 
and  53. 

'=•'7^  K-a-2Nv)p.]  (62) 

=0^  (53) 

where 

b^  =  M(N  +  l)v  +  [MN  + 1  -  (M  +  N)v}Yi. 
bj  =  M(N  +  l)v  -(MN  +1  -(M  +  N)vtY2, 

Cz  =  2^^[MN  +  l-(M  +  l)(N*l)vK  -(l-2Nv)pJ 

Expressions  for  stress  and  strain  in  the  elastic  zone  were  developed 

ar{R)  =  -2(l-Nv)Pb+au  (50) 


2Gee(R)  =  -2(N  -l)vpb-a. 


(51) 


In  Equation  53.1, 

,  -Cz-b,c,^  1  >|_  B.  ,  _b,c,-c,r  1  B, 

^  b, -b^  Ir^'-'J  R^‘-'’  ^  b,-b^  R-""-' 


(53.1) 


expressions  were  found  for  the  unknown  constants  Ai  and  A2  in  terms  of  the 
still  undetermined  radius  R  to  the  elastic-plastic  boundary. 
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Note  that  to  this  point  we  have  merely  reiterated  conditions  developed 
in  Section  3.  For  Case  II.  substitution  of  the  negative  of  the  internal 
pressure  -pa  for  Or  and  the  radial  and  circumferential  strains  from  Equations 
47  and  48  into  Equation  42  yields 


bjAja^'"*  +b2A2a'^'‘* 


.1-2NV 


(77) 


where  the  terms  bi,  b2.  and  C3  are  as  defined  in  Equations  53  and  55. 


Using  the  definitions  of  Ai  and  A2  from  Equation  53.1  gives  the 
following  equation  from  which  the  radius  R  to  the  elastic-plastic  boundary 
can  be  calculated. 


1-2NV 


(78) 


Note  that  the  corresponding  equation  in  Reference  4  (Equation  38)  appears  to 
contain  a  typographical  error;  the  second  term  exponent  should  be  1+a. 

As  long  as  Or  <  Oz  (the  outer  plastic  zone  jf  Figure  2),  radial  and 
circumferential  strains  can  be  computed  from  Equations  47  and  48.  Equations 
42  and  43  can  then  be  used  to  get  the  corresponding  stresses. 

4.2  INCREASED  LOADING. 

As  with  Case  I.  it  is  possible  to  increase  the  external  load  to  a 
magnitude  such  that  the  out-of-plane  stress  will  equal  the  radial  stress. 
This  condition  results  in  formation  of  the  inner  plastic  zone  of  Figure  2. 
In  this  situation,  the  radius  R  can  no  longer  be  obtained  from  Equation  78. 

As  Pb  is  increased,  the  inner  plastic  zone  begins  to  form  when  the  out- 
of-plane  stress  first  equals  the  radial  stress  at  the  edge  of  the  opening, 
i.  e.,  at  r  =  a.  The  far-field  pressure  at  this  instant  is  defined  as  Pb, 
and  the  corresponding  radius  to  the  elastic-plastic  boundary  is  defined  as 
R'.  The  development  of  expressions  for  these  terms  is  described  in  the  fol¬ 
lowing  paragraphs. 

The  expressions  for  radial  and  circumferential  stress  in  the  outer 
plastic  zone  were  presented  in  Equations  42  and  43.  Because  stress 
continuity  is  required  across  the  boundary  between  the  two  plastic  zones. 
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these  equations  hold  at  r=a  when  Pb~Pb.  Upon  invoicing  the  definition  of  C3 
from  Equation  55.  Equations  42  and  43  become 


CjO,  =  [MN  + 1  -  (M  +  N)v]2Ge,  +  M(N  +  l)v2Gee  +  (M  -l)vo„ 

CjCe  =  N(M  +  l)v2Ge,  +MN2Ge0-[l-(M  +  l)v]a„ 

The  out-of-plane  stress  from  the  yield  condition  of  Equation  1  is 


(79) 

(80) 


(81) 


Combining  Equations  80  and  81  yields 

030,  =  (M  +  l)v2Ge,  +  M2Gee  -  i[l  -  (M  +  l)v  -  C3 ]o„  (82) 

At  the  edge  of  the  opening,  i.  e.,  at  r  =  a,  Or  =  Oz  ~  "Pa.  and  strains 
er(a)  and  ee(a)  may  be  obtained  from  Equations  47  and  48.  When  these  substi¬ 
tutions  are  made.  Equations  79  and  82  become 


-C3P.  =[MN  +  1-(M+N)v] 


^[M(N+l)v] 


/  .  yi-i  /  3  \-Yi-i 

Y,c;B3(^-j  +A 

/gyi-i  f 

-"A 


+  (M-l)va 


(83) 


-C3P.  =[(M  +  l)v] 


YiCi'Bj 


r  /  a  ^ 

f  a ) 

-Ti-l 

+M 

f - 

w 

+  A 

(84) 


where  the  constants  83  and  84  are  obtained  by  rearranging  the  expressions  for 
81  and  82  (Equations  53.1) 


C.B3  = 


b.-b, 


Cad-Nv) 

(N-l)V 


-b. 


(85) 


qB, 


Cl  b  C3(1-Nv) 
b3-b,['  (N-l)v 


(86) 
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The  term  c'  is  defined  as  the  value  of  ci  (from  Equation  52)  when  Pb=Pi 


ci 


2(N  -l)v 
1-2NV 


-(l-2Nv)p;] 


(87) 


Equations  83  and  84  may  be  solved  simultaneously  and,  after  invoking 
the  definition  of  c' f rom  Equation  87,  we  obtain 


1 


^  {2vo„-a-2Nv)p,}b,-{o, 
B/{2vo„-a-2Nv)p.}b3-{o„ 


-(l-2Nv)p.}{M^(Mn)vYj 
-(1-2Nv)p.}{M-(M  +  1)vyJ 


-rrifi 


(88) 


and 


r 


b 


-  Ou 
1-2NV 


\l-2Nv 


- - P 

>11.. 


2(N-l)v 


b,B3 


(89) 


Since  p^  is  the  minimum  far-field  pressure  that  will  cause  a  second  plastic 
zone  to  form,  this  equation  can  be  used  to  determine  whether  one  or  two 
plastic  zones  will  form. 

The  maximum  radius  at  which  the  radial  and  out-of-plane  stresses  are 
equal  was  defined  earlier  as  R.  An  expression  for  the  ratio  of  R  to  R  can  be 
developed  in  a  fashion  similar  to  that  employed  in  the  development  of 
Equation  56  for  Case  I. 

The  stresses  at  r  =  R  can  be  expressed  by 

a,(R)  =  o,(R)  (90) 


As  stated  above,  the  out-of-plane  stress  from  the  yield  condition  of  Equation 
1  is 


Oz  = 


% 

N 


(81) 


Circumferential  stress  (at  r  =  R)  from  Equation  43  is  substituted  in  Equation 
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81,  which  is  then  equated  to  the  radial  stress  from  Equation  42  to  get 


[M(N  -  2  V)  -  Nv  +  a  -  v)]2Ge,  +  M[Nv  -  a  -  v)]2Ge9  =  M(1  -  2  v)o„  (91) 

Substitution  of  strains  from  Equations  47  and  48  in  Equation  91  gives 


b3AiR^‘''  +  b^A2R‘^*'‘ 


-2v 

-2Nv 


(92) 


Equation  92  is  identical  to  Equation  55  and  the  terms  bs.  b4  and  C3  were 
defined  in  conjunction  with  that  equation.  Using  the  definitions  of  Ai  and 
A2  from  Equations  53.1  yields  an  expression  for  the  desired  ratio 


—  +  b.B, 

U  J 

-rj-i 


l-2v 

=  c, - 

'1-2NV  “ 


(93) 


which  is  identical  to  Equation  56.  As  with  Case  I,  R  remains  unknown  pending 
consideration  of  the  region  where  Cr  =  Oz*  Also,  as  with  Case  I,  both  yield 
conditions  (Equations  57  and  58)  are  applicable  when  Or  =  Cfz  (the  inner 
plastic  zone  in  this  case). 


fi  =  Cq  -NOr  +au  =  0 

(57) 

f2  =  a0-Naz+au  =0 

(58) 

Using  Equation  57  and  the  equilibrium  expression  of  Equation  13 

r^^  +  ar-CTe  =  0  (13) 

dr 


Radial  stress  can  be  expressed  by  Equation  60. 


Or  =- 

Pa  ^ 

"  N-1 

(60) 


In  order  to  satisfy  radial  stress  continuity  at  r  =  R,  we  again  set  Equation 
60  equal  to  Equation  42  and  substitute  strains  from  Equations  47  and  48,  all 
quantities  being  evaluated  at  r  =  R.  As  with  Case  I,  this  yields  Equation 
62. 


biBi 


+  ^282 


=  C-, 


N(l-2v) 

(N-l)(l-2Nv) 
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Pa-^ 


N-1 


Oy  -C3 


(62) 


Since  R/R  was  determined  by  Equation  93,  R  can  be  computed  from  Equation  62, 
as  was  done  in  Case  I,  and  R  can  then  be  obtained  explicitly. 


The  flow  rule  expressed  by  Equations  64,  65,  and  66 


p(p)  _  iL  .  y  ^2  -  uT, 


;(p)  _  -V  afi  ,  afj 


p(p)  _  il  .^^  +  1  =  _m5l 


(64) 

(65) 

(66) 


again  leads  to  the  plastic  strain  expression  of  Equation  67: 
g(r )  =  E^f  ^  =  0 

At  r  =  R,  the  plastic  strain  expression  of  Equation  68  follows 
g{r)  = 

which,  with  Equations  40  (which  hold  in  the  outer  plastic  zone), 

+  MEq  '  =  0, 


e(p)=0 


leads  to  Equation  69; 


E).'”  +  MEg”’  +  e)'’’  =  0 


(67) 


(68) 


(40) 


(69) 


Invoking  the  plane  strain  constraint,  radial  strain  can  be  expressed  by 
Equation  70 


Er  =  -MEe  +MEe*’  +E^*’  +  e)*’ 


(70) 


Since  Or  =  Oz  in  the  inner  plastic  zone,  the  Case  I  derivation  of  Equations 
71  and  72  also  holds. 


Ee'*'  =Ee)*’  =  ll-(N  +  l)vla,  +VO, 


(71) 
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e4®^  =(N-2v)ar  -Ou 


(72) 


As  with  Case  I,  a  differential  equation  for  the  circumferential  stress 
in  terms  of  r  is  given  by  Equation  73  when  elastic  stresses  from  Equations  71 
and  72  and  radial  strain  from  Equation  70  are  substituted  in  the 
compatibility  expression  of  Equation  28. 


MN42-2(M>N-^1)v 
(N-l)(l  +  v) 


[(N-l)p. 


.N.H-l 


N-l 

d 


^(M4  2)(1-2v)c, 
(N-Dd-^v) 


(73) 


Again.  Equation  73  can  be  solved  by  direct  integration  from  r  to  R  to 
get  the  expression  for  circumferential  strain  in  Equation  74.  As  before, 
e0(R)  is  obtained  from  Equation  48.  Radial  strain  is  given  by  Equation  75 
which  was  obtained  from  Equations  70,  71.  and  72. 


2Gee  =  2Gee(R)- 


IT' 


HN-t-2-2(M-t-N-^l)v 
(M  +  N)(N-l)(l  +  v) 


[(N-l)p. 


+ 


(-T 


(M-^2)(l-2v)a„ 
(N-l)(M  +  l)(l  +  v) 


(74) 


2G£,  = -2GMe0  + 


MN  +  2-2(M  +  N+l)v 
- a, 


M-2v 
1  +  V 


(75) 


Florence  and  Schwer  (Reference  4)  determined  that  for  Case  II  with  an 
associated  flow  rule,  an  inner  plastic  zone,  where  dj  was  the  intermediate 
principal  stress,  did  not  form.  Since  our  equations  for  stress  in  the  inner 
plastic  zone  are  identical  to  theirs  (i.e.,  they  do  not  involve  M),  their  de¬ 
termination  holds  for  the  case  of  arbitrary  dilatancy. 
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SECTON  5 


APPLICATION  OF  CLOSED-FORM  SOLUTIONS 

In  Section  5.1,  procedures  for  applying  the  solutions  developed  in  the 
foregoing  sections  are  discussed.  Some  numerical  examples  are  presented  in 
Section  5.2. 

5.1  PROCEDURES. 

The  solutions  for  problems  of  axi symmetric  compression  of  Mohr-Coulomb 
materials  with  arbitrary  dilatancy  are  quite  complex,  as  are  the  procedures 
for  applying  them.  Consequently,  we  have  developed  a  "road  map"  that  is  in¬ 
tended  to  aid  the  reader,  first  in  determining  the  appropriate  solution  case 
(i.  e..  Case  I  or  Case  II  yielding),  and  then  in  selecting  the  proper  equa¬ 
tions  to  define  the  stress  and  strain  fields. 

The  proper  solution  case  may  be  determined  from  the  conditional  in¬ 
equalities  listed  earlier  in  Table  2.  These  inequalities  involve  only  the 
quantities  Pa,  Pb.  v,  Oy,  and  cp  (or  N).  These  and  the  other  parameters  that 
must  be  known  (either  given  or  assumed)  at  the  beginning  of  the  problem  are 
listed  in  Table  3.  The  process  is  detailed  in  Table  4.  The  decision  logic 
of  Table  4  is  illustrated  schematically  in  Figure  3.  Note  that  the  process 
of  determining  the  proper  yield  condition  is  independent  of  the  dilatancy 
parameter  M. 

As  may  be  seen  from  Table  4  and  Figure  3,  there  are  four  possible 
outcomes:  (1)  the  Mohr-Coulomb  medium  yields  everywhere.  (2)  the  medium  is 
elastic  everywhere,  (3)  the  medium  yields  under  Case  I  conditions,  or  (4) 
the  medium  yields  under  Case  II  conditions.  In  both  Cases  I  and  II,  the 
far-field  medium  remains  elastic.  Depending  upon  which  condition  prevails, 
the  reader  is  directed  to  Table  5,  6,  7,  or  8  for  the  appropriate  step-by- 
step  procedure.  Because  of  the  length  of  these  tables,  they  are  located  at 
the  end  of  this  section. 

As  indicated  in  Table  4  and  Figure  3.  if  under  Case  I,  (N+l)v  ^  1, 
only  a  single  plastic  zone  will  form,  regardless  of  the  value  of  Pb^Pb.  This 
condition  follows  from  the  fact  that  Oz  will  be  the  intermediate  principal 
stress  throughout  the  plastic  zone  when  (N+l)v  ^  1.  On  the  other  hand,  if 
(N+l)v  <  1,  the  three  plastic  zones  shown  in  Figure  1  will  form  as  Pb  is 
increased. 
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Table  3.  Known  parameters. 


The  following  parameters  must  be  either  given  or  assumed: 

Unconfined  compressive  strength  Oy 

r  .  . .  ,  ,  l  +  sin(p . 

Friction  angle  (p  (or  N  = - —) 

l-sin(p 

Poisson's  ratio  v 

Arbitrary  dilatancy  factor  M,  1<M<N 
Internal  pressure  Pa 
External  pressure  Pb 
Shear  modulus  G  (or  E  =  2G(l+v)) 

Interior  radius  a  (not  required  when  results  are  normalized  with 
respect  to  radius ) 
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Table  4.  Process  for  determining  proper  step-by-step 
procedure  to  follow. 

If  Nv  <  ^  and  Pa  >  ^  .  then  the  far-field  yields  everywhere  and  the 

solution  is  not  valid--quit. 

Otherwise,  coirpute  (N  +  l)v 

If  (N+l)v^l,  Case  I  yielding  will  occur  for  Pb^Pb*  Use  the  single-plastic- 
zone  procedure  of  Table  6. 

Otherwise  (i.e.,  (N+l)v<l),  compute 

If  Pa  <  ^  yielding  will  occur  for  Pb^^Pb-  Use  the  three- 

plastic-zone  procedure  of  Table  7. 

Otherwise,  Case  II  yielding  will  occur  for  Pb^^Pb.  Use  the  procedure  of 
Table  8. 
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In  the  three-plastic-zone  case,  only  a  single  plastic  zone  will  form 
when  the  far-field  pressure  is  small,  i.  e.,  when  p^<Pi,<Pb.  However, there 
will  not  be  a  situation  under  Case  I  where  two  plastic  zones  form.  This  is 
true  because  the  middle  and  outer  plastic  zones  begin  to  form  simultaneously 
when  pb  =  Pb  because  R  =  R  =  R  at  that  instant. 

Under  Case  II  yielding,  the  general  situation  is  two  plastic  zones  as 
shown  in  Figure  2.  However,  as  with  Case  I,  only  a  single  plastic  zone  will 
exist  at  low  far-field  pressures,  i.  e.,  when  Pb^Pb<Pi. 

In  the  step-by-step  procedures  (Tables  5  through  8)  the  equations  de¬ 
veloped  in  the  foregoing  sections  are  arranged  in  order  to  allow  one  to  ex¬ 
peditiously  determine  the  stress  and  strain  fields  in  the  medium  surrounding 
the  tunnel.  Appropriate  amplifying  comments  have  also  been  included  in  the 
tables.  In  general,  each  procedure  begins  with  the  computation  of  required 
constants,  followed  by  the  determination  of  the  plastic-plastic  and  elastic- 
plastic  radii,  and  concludes  with  the  equations  necessary  to  completely 
define  the  stress  and  strain  fields  in  each  zone.  For  the  reader  who  is  in¬ 
terested  in  determining  only  the  circumferential  strain  at  the  tunnel  wall 
(equivalent  to  the  negative  of  tunnel  closure  (AD/D)  in  axi symmetric 
problems  of  this  type),  the  special  case  for  eeat  r  =  a  is  presented  at  the 
end  of  each  table. 

Equation  numbers  from  the  foregoing  sections  have  been  repeated  in 
Tables  5  through  8.  A  primed  equation  number  indicates  that  the  equation 
has  been  modified  in  some  way.  Where  two  equations  have  been  combined,  the 
numbers  of  both  are  given.  Equations  that  appear  in  the  tables  without 
numbers  are  unique  to  the  step-by-step  procedures.  In  the  tables,  we  have 
continued  the  use  of  twice  the  shear  modulus  times  the  strain,  rather  than 
expressing  the  strains  explicitly. 

5.2  NUMERICAL  EXAMPLES. 

Two  numerical  examples  are  presented,  using  the  input  parameters  spec¬ 
ified  by  Florence  and  Schwer  in  References  3  and  4.  In  addition  to  indicat¬ 
ing  some  of  the  effects  of  arbitrary  dilatancy,  these  exanples  illustrate 
the  process  for  determining  the  appropriate  yield  case. 

5.2.1  Case  I  Example. 

For  this  example,  input  parameters,  taken  from  Florence  and  Schwer 
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(Reference  3),  are: 

Unconfined  compressive  strength,  Ou  =  2,000  psi  (13.8  MPa) 

Friction  angle,  (p  =  30®,  N  =  3 

Poisson's  ratio,  v  =  0.2 

Internal  pressure,  Pa  =  500  psi  (3.45  MPa) 

Far-field  pressure,  Pb  =  6,000  psi  (41.4  MPa) 

Shear  modulus,  G  =  10®  psi  (6,895  MPa) 

Interior  radius,  a  (not  specified— results  presented  as  functions  of 
r/a ) 

Three  different  values  of  the  arbitrary  dilatancy  parameter  are  used.  These 
are  (1)  M  =  N,  (2)  M  =  1.  and  (3)  M  =  (N  +  l)/2. 

The  first  step  is  to  determine  the  appropriate  step-by-step  procedure. 
Starting  at  the  top  of  Table  4,  Nv  =  0.6  >  so  the  solution  procedure  is 
valid. 


Next,  (N  +  l)v  =  0.8.  Therefore,  (N  +  l)v  <  1.  In  addition. 


1  -(n  +  l)v 


=  800 


and 


p,  =  500<  800 

which  indicates  that  the  Case  I  yield  condition  applies,  so  the  procedures 
of  Table  7  should  be  used. 

The  first  step  in  Table  7  is  to  determine  whether  the  specified  value 
of  Pb  will  cause  yielding  to  occur  and,  if  so.  whether  one  or  three  plastic 

zones  will  form.  Therefore, 

p,  =  6,000  >i[(Nn)p.+aJ  =  2,000 

which  indicates  that  yielding  occurs.  Next 

- — ^ - -  =  5,000  <  p  =  6,000 

2[l-(N  +  l)v]  " 

Therefore,  Pb  is  sufficiently  large  to  cause  the  formation  of  three  plastic 
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zones. 


For  the  upper-bound  case  where  M  =  N,  calculated  radius  ratios  are  (1) 
R/a  =  V2,  (2)  R /a  =  1.527,  and  (3)  R/a=1.598.  Stresses  and  strains  for  this 
case  are  plotted  in  Figure  4.  Except  for  differences  in  scale,  these  plots 
are  identical  to  the  Florence  and  Schwer  results  (Figures  2  and  3  of 
Reference  3).  Contrary  to  the  established  sign  convention,  we  also  plotted 
all  compressive  stresses  in  Figure  4  (and  subsequent  figures)  as  positive  to 
facilitate  comparisons  with  the  Florence  and  Schwer  results.  Tensile 
strains  are  shown  as  positive,  both  here  and  in  Reference  3. 

For  the  lower-bound  case  where  M  =  1,  the  calculated  radius  ratios  are 
(1)  R/a  =  V2.  (2)  R /a  =  1.523,  and  (3)  R/ a  =1.596.  For  the  intermediate  case 
where  M  =  (N  +  1  )/2  =  2,  the  ratios  are  (1)  R/a  =  V2,  (2)  R/a=1.526,  and  (3) 
R/a=1.598. 


The  ratio  R/a  remains  constant  as  M  is  varied  since,  as  may  be  seen  in 
Equation  38,  it  is  not  a  function  of  M.  The  other  two  ratios  change 
slightly  when  M  is  increased  from  M  =  1  to  M  =  N. 

Stresses  and  strains  for  all  three  cases  are  plotted  in  Figure  5. 
Note  that  the  horizontal  scale  is  more  than  double  that  of  Figure  4.  For 
this  example,  there  are  no  discernible  differences  in  the  stress  plots. 
Strains  differ  significantly  at  the  edge  of  the  opening,  but  are  nearly 
identical,  regardless  of  the  value  of  M,  beyond  r/a  of  about  1.5. 

5.2.2  Case  II  Example. 

Input  parameters  for  this  example  were  taken  from  Florence  and  Schwer 
(Reference  4): 

Unconfined  compressive  strength,  Ou  =  3,000  psi  (20.7  MPa) 

Friction  angle,  cp  =  19.5°,  N  =  2 

Poisson's  ratio,  v  =  0.2 

Internal  pressure,  Pa  =  2,000  psi  (13.8  MPa) 

Far-field  pressure,  Pb  =  6,000  psi  (41.4  MPa) 

Shear  modulus,  G  =  0.5  x  10®  psi  (3,447  MPa) 

Interior  radius,  a  (not  specified--results  presented  as  functions  of 
r/a ) 

Again,  three  different  values  of  the  arbitrary  dilatancy  parameter  are  used: 
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STRAIN  (PERCENT)  STRESS  (KSI) 

0.000  1.500  3.000  4.00  8.00  12.00 


Figure  5.  Computed  stresses  and  strains  for  Case  I  example  with  different 
values  of  M. 
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(1)  M  =  N.  (2)  M  =  1.  and  M  =  (N  +  l)/2. 


Beginning  at  the  top  of  Table  4,  Nv  =  0.4<|^.  However, 


—  =  15, 000  >p  =2,000 
(1-2NV) 


so  the  solution  procedure  is  valid.  The  next  condition  is  (N+l)v  =  0.6. 
Therefore,  (N  +  l)v<l.  In  addition. 


ygu 

l-(N  +  l)v 


1,500 


and  p,  =  2,000  >1,500,  which  indicates  that  the  Case  II  yield  condition 
applies,  so  the  procedures  of  Table  8  should  be  used. 

The  first  step  in  Table  8  is  to  determine  whether  the  specified  value 
of  Pb  will  cause  yielding  to  occur.  Therefore, 

p  =6,000  >  =4,167 

2(1 -Nv) 


which  indicates  that  yielding  occurs. 

Next  it  is  necessary  to  determine  whether  one  or  two  plastic  zones 
will  form. 


In  the  case  of  M  =  N, 
p;  =  4,529<p,  =  6,000 

so  the  two-plastic-zone  procedure  applies. 

For  the  upper-bound  case  where  M  =  N,  calculated  radius  ratios  are  (1) 
R'/a  =  1.100,  (2)  R/a  =  1.190,  and  (3)  R/a  =  1.555.  Stresses  and  strains 
for  this  case  are  plotted  in  Figure  6.  Except  for  differences  of  scale, 
these  plots  are  identical  to  the  Florence  and  Schwer  results  (Figures  2  and 
3  of  Reference  4).  Florence  and  Schwer  state  that  they  used  a  shear  modulus 
G  =  10^  psi  (6,895  MPa).  However,  their  plotted  strains  are  consistent  with 
2G  =  10^  psi  (6,895  MPa). 

The  two-plastic-zone  procedure  also  applies  for  the  other  two  values 
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STRAIN  (PERCEMT)  STRESS  (KSI) 

10  0.000  0.500  3.00  6.00  9.00 


Figure  6.  Computed  stresses  and  strains  for  Case  II  example  with  M  =  N. 
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of  M.  For  the  lower  bound  case  where  M  =  1,  the  calculated  radius  ratios 
are  (1)  RVa  =  1.114,  (2)  R/a  =  1.175,  and  (3)  R/a  =  1.566.  For  the 
intermediate  case  where  M  =  (N+l)/2  =  1.5,  the  ratios  are  (1)  R'/a  =  1.104, 
(2)  R/a  =  1.185,  and  (3)  R/a  =  1.558.  Changes  in  the  radius  ratios  are 
quite  small.  R'/a  and  R/a  decrease  by  1.3  and  0.7  percent  respectively 
while  R/a  increases  by  1.3  percent  as  M  is  increased  from  M  =  1  to  M  =  N. 

Stresses  and  strains  for  all  three  cases  are  plotted  in  Figure  7. 
Again,  note  that  the  horizontal  scale  is  more  than  double  that  of  Figure  6. 
Stresses  vary  only  slightly  with  changes  in  M.  Radial  strains  differ 
significantly  at  the  tunnel  wall,  but  circumferential  strains  differ  only 
slightly.  Both  stresses  and  strains  are  nearly  identical,  regardless  of  the 
value  of  M,  beyond  r/a  of  about  1.5. 
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Table  5.  Procedure  for  elastic  condi c ;  . 


The  strain  and  stress  fields  are  given  by 

2Ge,  =-(l-2v)p,+(p, -p.)-^  (6) 

2GEe  =  -(l-2v)p,  -(p,  -pj-^  (6) 

r 

gZ 

=-Pb+(Pb-P.)-T 

r 

gZ 

=  -Pb-(Pb -P.)— 

r 

a,  = -2vp,  (5) 

TUNNEL  CLOSURE 

Y  =  -e9(a)  =  “[P.-2Pb(l-v)l  (6') 
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Table  6.  Procedure  for  Case  I  yielding— single  plastic  zone. 


For  low  values  of  Pb  ( Pb  ^  2[(N+l)p, +cru] ).  use  the  procedure  for  elastic 
conditions  of  Table  5.  Otherwise: 


Compute  R 

2  (N-l)pb-»Ou 

UJ  N+1  {N-l)pj  +  au 


(34) 


then  R  =  af-L 

yaj 

Plastic  Zone  a <  r  < R 


N-1 


N-l 


—  +■ 


N-l 


2Ge^*'=[N-(N  +  l)v]o,-(l-v)a, 


H-.X 


2Gee  =  2Ge^*>  +  2Ge^^' 

2Ge,  =  [1  -  (N  +  l)v]a,  +  va„  -  M2Ge^'’^ 


Oe  =  Na,  -cj, 

o,  =  v(o,  +Oe) 


(14) 

(26) 

(30') 

(25.17) 

(7') 

(24) 


Pa  is  required  in  the  elastic  zone.  Therefore, 

g^(R)  =  JPbl^  (33) 

N+1 

Pi  =  -<^r(R) 
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Table  6.  Procedure  for  Case  I  yielding— single  plastic  zone 
(Concluded). 


Elastic  Zone  r  >  R 


2G£  =-(1-2v)p,+(P,-p:)^ 

(6') 

r 

r2 

2Gee  =  -(l-2v)pb-(Pb-p*,)— 

r 

(6’) 

tJr  =-Pb+(Pb-Pl)^ 

(5’) 

r 

r2 

O0  =  -Pb-(Pb-Pl)-T 

(5') 

r 

CT,  =  -2vp, 

(5) 

TUNNEL  CLOSURE 


D 


-ee(a)  = -^[{N -(N +l)v}p. -d  -  v)aj 

cva 


2G 


N  +1 


.R 


(26,30') 
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Table  7.  Procedure  for  Case  I  yielding— three  plastic  zones. 

For  low  values  of  Pb  ( Pb  ^ ^^e  the  procedure  for  elastic 
conditions  of  Table  5.  Next,  determine  whether  there  will  be  one  plastic  zone 
or  three.  If  Pb<  (N^l)v]*  ^  single  plastic  zone  will  form  and  the 

procedure  of  Table  6  may  be  used.  Otherwise: 

Compute  required  constants 


MN  +  1-(M  +  N)v 

(45) 

,„.(N-M)v  ,  l(N-M)'v=^  1 

2MN  ^  \  (MN)^ 

(49) 

''  2MN  (MN)‘  P' 

(49) 

1-2NV  “ 

(49) 

C  [o^ -(l-2Nv)p,] 

(52) 

bj  =M(N  +  1)v  +  {MN  +  1-(M  +  N)vIy, 

(53) 

b^  =  M(N+1)v-(MN  +  1-(M  +  N)v]y2 

(53) 

c^  =  + 1  - (M  +  1)(N  +  l)v][a„  - (1  -  2Nv)p J 

(53) 

g  c^-b^c,^ 

bj-b^  Pi’^z 

(53.1) 

bj  =  M[(N  +l)v  -1]  +  [MN  +  1  -(2M  +  N  +1)v|Yi 

(55) 

b,  =  Ml(N  +l)v  -1]  -[MN  +  1  -(2M  +  N  ^DvIy^ 

(55) 

C3  =[MN  +  l-(M+l)(N  +  l)v] 

(55) 
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Table  7.  Procedure  for  Case  I  yielding-three  plastic  zones 
(Continued). 


Now  determine  radii  R,  R,  R 

_i 

(i-2v)qu _ 

{l-(N+l)v}{(N-l)p.+c„} 


(38') 


Then 


R  =  a 


The  ratio  —  may  be  determined  by  iteration. 
R 

solve 


Begin  by  assuming  —  =  0.9  and 


R 

R 


l-2v  .01” 
'1-2NV  “  '  Hr 


b,B3 


-Y,-l 


(56') 


p 

Compute  —  from 
a 


\^J 


rV.-‘ 


r; 


r 


+  c. 


N(l-2v) 

(N-l)(l-2Nv)' 


C3  P.+ 


N-1 


(62’) 


then 


R  =  d 


-nN-1 


.a. 


H-l 


and 


v'' ; 


Compute 


(63) 
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Table  7.  Procedure  for  Case  I  yielding-three  plastic  zones 
(Continued). 

pi  is  required  in  the  elastic  zone.  Therefore  in  the  outer  plastic  zone, 
compute  strains  and  stresses  at  r  =  R  (the  outer  boundary  of  the  zone) 


2Ge,(R)  =  Yl^R^'■'  -YzAzR'^''^  +  A 

(47') 

2Gee(R)  =  AjR^‘'‘  +  A2R'^*‘‘  +  A 

(48') 

a,  (R)  =  — [{MN  + 1  -  (M  +  N)v}2Ge,  (R)  +  M(N  +  l)v2Gee(R)  +  (M  -l)va„  ] 

^3 

(42') 

p;=-o,(R) 

Inner  plastic  zone  strains  depend  on  circumferential  strains  at  r  =  R  and  r=R 
so  those  are  computed  next. 

2Gee(R)  =  +  A  (48’) 


2Gee(R)  =  2Gee(R). 


^^Y*‘^MN  +  2-2(M  +  N  +  1)v 


V' 


(M  +  N)(N-l)a+v) 


[(N-l)p.+aJ 


u 


-1 


--sH-l 


(M^2)(l-2v)q„ 


(N-l)(M+l)(l  +  v) 


vR/ 


(74-) 


The  plastic  component  of  ee(R)  is  also  required  in  the  inner  plastic  zone.  It 
may  be  determined  as  follows 


l-(N  +  l)v 


(35) 


2Ge‘*’{R)  =  [N  -(N  +1)v1o,(r)  -(1  -  v)a„  (26') 

2Ge''”(R)  =  2G£e(R)-2Ge‘*>(R)  (19’) 


With  the  information  developed  as  described  above,  stresses  and  strains  at  any 
point  may  be  readily  determined.  For  a  point  at  any  radius  r  >  a,  compare  r 
to  R,  R,  and  R  to  determine  its  location;  then  proceed  to  use  the  appropri¬ 
ate  formulas  below. 
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Table  7.  Procedure  for  Case  I  yielding— three  plastic  zones 
( Continued). 


Inner  Plastic  Zone 


Or  =  -  [p, 

L  N-1 


(a  ^  r  <  R) 


06  =  No,  -a„ 


o,  =  v(a,  +ae) 

2Ge‘*^  =  [1  -  (N  +  l)v]o,  +  vOu 
2Ge'*^=[N-(N+l)v]a,-a-v)o„ 


/  ■ 

2Ge<^^=2Ge'^'(R)  - 


(N>l)(l-v) 
M  +  N 


[(N-Dp.+oJ  - 


-xN-l 

R' 


2Ge6  =  2Ge^*^  +  2Ge^'’^ 

2Gei‘’'  =  -M2Ge‘6^’ 

2Ge,  =  2Get*^  +  2Get'’^ 

Middle  Plastic  Zone  (R<r<R) 


2GEe  =  2Ge6(R)[- 


MN  +  2-2(M  +  N  +  l)v^^., 

+ - ^  (N  -l)p,  +  o J 

(M  +  N)(N-l)(l  +  v)  * 


Vr  J 


-1 


N-l 


(M  +  2)(l-2v)o„ 


(N-l)(M  +  l)(l  +  v) 


(3"- 


(14) 

(7’) 

(24) 

(25) 

(26) 

(76) 

a9‘) 

(16’) 

(19') 


(74) 


o.  uor  MN  +  2-2(M  +  N+l)v  M-2v 

2Ge,  =  -M2G0e  + - - —o, - o, 

'  ®  1+v  '  1+v 

N-l  _ 


P.  "^TTr 

- 

*  N-1 

vay 

N-1 


Oe  =  Na,  -o, 
o,  =o. 


(75) 

(14) 

(7') 
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Table  7.  Procedure  for  Case  I  yie1ding--three  plastic  zones 
(Continued). 


Outer  Plastic  Zone  (R<r<R) 

2Ge,  =  '  "YzAzr'^*"'  +  A 

2Gee  =  Ajr^*''  +  A^r'^'’^  +  A 

a,  =  — [{MN  +  l-(M  +  N)v}2Ge,  +M(N  +  l)v2Gee  +  (M-l)vaJ 

^3 

Oe  =  — [N(M  +  l)v2Ge,  +  MN2Gee  -  {1  -  (M  +  l)v}a„  ] 

^3 

Elastic  Zone  (r>R) 

2Ge,  =-(l-2v)p,+(p,-p;)^ 

2Gee  =  -(l-2v)p,-(p,-p;)^ 

tJr  =-Pb-^(Pb-Pl)^ 
r 

r2 

06=-Pb-(Pb-pi)^ 
a,  =  -2vp, 


(47) 

(48) 

(42’) 

(43') 

(!') 


(6‘) 

(6') 

(S') 

(S') 

(S) 
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Table  7.  Procedure  for  Case  I  yielding— three  plastic  zones. 
(Concluded). 


TUNNEL  CLOSURE 


Compute  constants  P*.  Yi.  Y2.  A.  bj.  b2.  b3,  b4,  c^,  C2.  C3,  Bi,  and  B2  as 
shown  at  the  beginning  of  the  table.  Then  coirpute  R,  R,  and  R,  followed  by 
and  A2. 


Now  confute  tangential  strain  components  as  follows. 
2Gee(R)  =  AiR^‘‘'  +  A2r^‘''+A 


fn\  or  /ox  ( R  MN  + 
2Ge9(R)=2Gee(R)-|^|-J  +-^ 


MN*2-2(M-»N-^l)v 
N)(N-l)(l  +  v) 


[(N-l)p.+aJ 


ir^r  ,1 

h- 

lUJ  J 

(M  +  2)a-2v)a 


(N-l)(M  +  l)(l  +  v) 

2Ge'*\R)  =  -(i.v)a, 

®  l-(N  +  l)v  “  “ 

2Ge|,'’'(R)  =  2Ge9(R)-2Ge‘*>(R) 
2Ge‘*'(a)=-  [N -(N  +  l)v]p, -a-v)o„ 


“ 

\ 

N«1 

R 

-1 

/-nN-1 
I  D  \ 


2Ge;'”(a)  =  2Ge^'”(R) 


(N^l)(l-v) 
M  +  N 


[(N-l)p.  +  aJ 


/■  - 


/  -  xN*! 


Then 


D 


-ee(a)  =  ~[2Ge^*'(a)  +  2Ge''»(a)] 


(AS-) 


(74') 

(26') 

(19’) 

(26') 

(76') 

(19') 
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Table  8.  Procedure  for  Case  II  yielding. 


p  +  c 

For  low  values  of  Pw  (p.  <— ^ — use  the  procedure  for  elastic 

2(1 -Nv) 

conditions  of  Table  5. 

Determination  of  the  number  of  plastic  zones  is  not  as  straight  forward 
as  in  Case  I  yielding  because  calculation  of  p^  from  Equation  89  is  required. 
It  is  recommended  that  all  of  the  following  constants  be  computed  first  since 
most  of  them  are  required  in  the  calculation  of  p^ 


MN  +  1-(M  +  N)v 


'  ^  2MN  ^ 


j(N-M)^v^ 

1 

(  (MN)^ 

(M-N)^v^ 

1 

(MN)^ 

*p' 

^  ^  2MN  ^  ^  (MN)' 

A=-»^o. 

1-2NV  • 

K-a-2Nv)p.] 

bj  =  M(N  + 1)  V  +  [MN  + 1  -  (M  +  N)v1Yj 
b^  =  M(N  +  1)v-[MN  +  1-(M  +  N)v}y2 
bj  =  M[(N  +l)v  -1]  +  [MN  +1  -(2M  +  N  +1)v]y, 
b,  =  M[(N  +  l)v  -1]  -  [MN  + 1  -  (2M  +  N  +  Dvjy^ 
C3  =  [MN  +  l-(M-*-l)(N+l)v] 


c,  =2 


,  1-Nv 
■l-2Nv 


[MN  + 1  -  (M  +  1)(N  + 1) v][o„  -  (1  -  2Nv)p,  ] 


B  = 

bj-b, 


B  =bi5^ 

^  b^-b^ 


(53) 

(53.1) 
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Table  8.  Procedure  for  Case  II  yielding. 
( Continued). 


83=^  (85’) 

c, 

8,=-^  (86') 

Cl 


Now  determine  whether  one  or  two  plastic  zones  will  form  by  computing  the 
ratio  a/R‘  and  then  P^ 


a 

F 


{2vq„  -  a-2Nv)p,}b3-{q„ 
{2vq„  -a-2Nv)p,}b2-{q„ 


-a-2Nv)p,}{M^(Mn)vYj 

-a-2Nv)p.}{M-(M+l)vYj 


1 

-Ti-ITt 


(88) 


and 


(89) 


If  Pb^P^,  use  the  single-plastic-zone  procedure  immediately  below. 
Otherwise,  proceed  directly  to  the  section  below  entitled  TWO  PLASTIC  ZONES. 

SINGLE  PLASTIC  ZONE 


The  ratio  a/R  may  be  determined  by  iteration.  Begin  by  assuming  a/R  = 
0.9  and  solve 


a 

R 


1-2NV 


1 

-Yi-l 


Then  R  =  a/(a/R).  Compute  constants 


(78') 


(63) 
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Table  8.  Procedure  for  Case  II  yielding 
(Continued). 

With  the  information  developed  as  described  above,  stresses  and  strains  at  any 
point  may  be  readily  determined.  For  a  point  at  any  radius  r  >  a,  compare  r 
to  R  to  determine  its  location;  then  proceed  to  use  the  appropriate 


formulas  below. 

Plastic  Zone  (a<r  <R) 

2Ge,  =yiA,r^‘‘'  +A 

(47) 

ZGEe  =  +  Ajr  +  A 

(48) 

o,  =  — [{MN  + 1  -  (M  +  N)v}2Ge,  +  M(N  +  l)v2Gee  +  (M  -l)vo,  ] 

^3 

(42’) 

Oe  =  — [N(M  +  l)v2Ge,  +  MNZGCe  -  {l  -  (M  +  l)vK  ] 

^3 

(43') 

(!') 

Elastic  Zone  (  r  >  R 1 

First  find  P4 

a,(R)  =  -2(l-Nv)p,+a„ 

(50) 

P*.  =  -o,(R) 

2Ge,  =-(l-2v)p,+(p,-p;)^ 

(6') 

2GEe  =  -(l-2v)p,-(p,-p;)^ 

(6') 

cy,  =-p,  +  (P5-p:)-^ 

(5‘) 

rZ 

09  =  -Pb-(Pb-Pl)  — 
r 

(5’) 

a,  =  -2vp, 

(5) 
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Table  8.  Procedure  for  Case  II  yielding 
(Continued). 

The  tunnel  closure  equation  is  provided  at  the  end  of  this  table. 
TWO  PLASTIC  ZONES 


Determine  the  radii  R  and  R.  The  ratio  R/R  may  be  determined  by  iteration. 
Begin  by  assuming  R/R  =  0.9  and  solve 


R 

R 


l-2v 

c, - a.  -b,B, 

'l-2Nv  “  '  ' 


b,B, 


1 

-Ti-l 


(56’) 


compute  R/a  from 


u  n  u  n  fR  V  N(l-2v) 

-b.B.  -  -b,B,  -  +c, - ^ — a. 

U  J  U  )  _ '  (N-l)(l-2Nv) 

3l. 

N-1 


Then 


C3  P.+ 


R  =  d 


(iT 


I 

N-l 


and  R  =  7= 


(R/R) 


Compute  constants 


^  RTT.-l  ’ 


-7,-1 


(62’) 


(63) 


P,  is  required  in  the  elastic  zone.  Therefore  in  the  outer  plastic  zone, 
compute  strains  and  stresses  at  r  =  R  (the  outer  boundary  of  the  zone) 


2Ge,(R)  =  Y^^R’'■■'  +  A 

(47’) 

2Ge,(R)  =  AiR'^'‘'  +  AjR''^''‘+A 

(48’) 
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Table  8.  Procedure  for  Case  II  yielding 
(Continued). 

Or  (R)  =  — [{MN  + 1  -  (M  +  N)v}2Ge,  (R)  +  M(N  +  l)v2Gee(R)  +  (M  -l)va„  ] 

^3 

P*.  =  -o,(R) 

Inner  Plastic  Zone  (r<R) 

2Gee(R)  =  AiR^'''+A2R’^*‘‘+A 


R 


2GEe  =  2Gee(R)-  - 


(M  +  N)(N-l)(l  +  v) 
(M  +  2)(l-2v)a„ 


N-l 


(N-l)(M+l)(l  +  v) 


Nrl 

-1 

2Ge,  = -M2Gee  + 


MN  +  2-2(M  +  N+l)v  M-2v 


1  +  v 


a.  -- 
"  1+v 


Or 


jEu-Ti 

N-lJU 


N-l 


+ 


Ou 

N-l 


Oe  =  Na,  -a„ 
o,  =  o. 

Outer  Plastic  Zone  (R<r<R) 

2Ge^  =Y,A,r^‘'*  "YA'" 

2Gee  =  Ajr^''^  +  A^r'^'"'  +  A 


(42‘) 


(48') 


(74) 


(75) 

(14) 

(7') 


(47) 

(48) 


Or  =  — [{MN  +  l-(M+N)v}2Ge,  +M(N  +  l)v2Gee+(M-l)vaJ 

^3 

Oe  =  — [N(M  +  l)v2Ge,  +  MN2Gee  -  {1  -  (M  +  l)v}a„  ] 

^3 


(42') 

(43') 

(!’) 
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Table  8.  Procedure  for  Case  II  yielding 
(Continued). 


Elastic  Zone  (r>R) 


2Ge,  =  -a-2v)p,+(p,-p;)\ 

r 

(6') 

2Gee  =  -a-2v)p,-(p,-p;)^ 

r 

(6') 

r2 

<yr  =-Pb+(Pb-Pl)-T 

r 

(5') 

r2 

<Je  =  -Pb-(Pb-Pl)— 
r 

(5') 

a,  =  -2vp, 

(5) 

TUNNEL  CLOSURE 

Compute  constants  Yi,  Yz,  A,  bj,  bj,  b3,  b4,  c^,  c^,  C3,  Bj,  B2,  B3, 
and  B4  as  shown  at  the  beginning  of  the  table.  Equations  88  and  89  can  then 
be  used  to  determine  if  one  plastic  zone  or  two  will  form. 

For  the  single-plastic-zone  case,  a/R  is  computed  from  Equation  78'  and 
constants  and  found  from  Equations  63.  Tunnel  closure  can  then  be 

determined  from 

When  two  plastic  zones  form,  R/R  and  R/a  are  computed  from  Equations 
56'  and  62',  respectively.  Then 

.  n  R 

and  R=7=— 

(R/ 


R  =  d  - 
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Table  8.  Procedure  for  Case  II  yielding 
(Concluded). 


Constants  and  A2  can  then  be  computed  from  Equations  63.  Then 


(«) 


and 


AD 


=  -e6(a)  =  -ee(R) 


MN  +  2-2(M+N  +  l)v 


2G(M  +  N)(N-l)(l+v) 


[(N-l)p.+oJ  -  -1 


(M  +  2)(l-2v)a 


2G(N-l)(M+l)(l  +  v) 


(74') 
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APPENDIX  A 

FLORENCE  AND  SCHUER  CASE  I  SOLUTION 


The  following  paper  (Reference  3)  was  published  in  the  International 
Journal  for  Numerical  and  Analytical  Methods  in  Geomechanics  by  John  Wiley  & 
Sons,  Ltd.  With  their  permission,  it  is  reproduced  here  in  its  entirety. 
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AXISYMMETRIC  COMPRESSION  OF  A  MOHR- 
COULOMB  MEDIUM  AROUND  A  CIRCULAR  HOLE 


A.  L.  FLORENCE  AND  L.  E.  SCHWER 
Poulter  Laboratory,  SRI  International,  Menlo  Park,  California  94025,  U.S.A. 


SUMMARY 

An  analytical  solution  is  presented  for  the  stress  and  strain  fields  in  a  Mohr-Coulomb  material  in  plane 
strain  around  a  circular  hole  when  it  is  compressed  by  an  axisymmetric  far-field  pressure.  It  is  shown  that 
several  solutions  arise  involving  one  to  three  plastic  zones  depending  on  the  values  of  Poisson's  ratio  and 
the  friction  angle.  The  solution  chosen  for  presentation  was  obtained  and  used  to  validate  the  functioning 
of  the  Mohr-Coulomb  yield  condition  that  was  added  to  the  NONSAP  finite  element  code.  Stress  and 
strain  field  comparisons  are  made. 


INTRODUCTION 

The  purpose  of  this  paper  is  to  provide  an  analytical  solution  of  a  rock  mechanics  problem  that 
can  be  used  to  validate  the  predictions  of  finite  element  codes.  The  problem,  which  is  of 
interest  in  itself,  is  the  determination  of  the  stress  and  strain  fields  surrounding  a  tunnel  in  a 
Mohr-Coulomb  rock  subjected  far  from  the  tunnel  to  an  axisymmetric  pressure.  Plane  strain 
conditions  prevail  and  an  internal  pressure  is  allowed  to  develop  at  the  hole  boundary.  Such  a 
pressure  may  be  considered  as  arising  from  a  metal  liner  at  yield  or  from  a  crushable 
back-packing  material  between  the  liner  and  the  rock.  To  simplify  the  analysis  and  to  provide 
a  defined  loading  history  this  internal  pressure  is  assumed  to  be  equal  to  the  far-field  pressure 
as  the  far-field  pressure  is  increased  (loading)  until  the  internal  pressure  reaches  its  final  value, 
whereupon  it  remains  constant  with  subsequent  loading. 

The  rock  behaviour  is  elastic-perfectly  plastic,  obeying  the  Mohr-Coulomb  yield  condition 
and  associated  flow  rule;  hence,  dilation  is  included.  It  will  be  seen  that  several  cases  can  arise 
that  depend  on  the  relative  magnitudes  of  the  principal  stresses,  and  these  relative  magnitudes 
depend  on  the  values  of  Poisson's  ratio  and  the  friction  angle.  One  of  these  cases  is  treated 
here,  chosen  to  involve  two  adjacent  faces  and  the  common  edge  of  the  yield  pyramid,  instead 
of  the  possible  case  involving  only  one  face,  to  provide  a  more  severe  test  for  the  code.  Three 
plastic  zones  arise,  shown  in  Figure  1,  that  correspond  to  the  two  faces  and  common  edge.  In 
this  work  the  finite  element  code  employed  is  NONSAP'  to  which  the  Mohr-Coulomb  yield 
condition  presented  in  Reference  2  has  been  added.  Comparison  of  analytical  and  code  results 
is  presented. 

The  motivation  for  the  code  validation  analysis  was  to  provide  support  for  the  experimental 
and  theoretical  investigations  of  Kennedy  and  Lindbergh  and  Senseny  and  Lindberg*  into  the 
response  of  buried  cylindrical  structures  to  quasi-static  loading. 

Initial  loading — outer  elastic  zone 

Let  p,  be  the  pressure  acting  on  the  hole  boundary  of  radius  a  and  let  Pb  be  the  pressure  at 
infinity  (far-field).  If  p,  =  pb  on  initial  loading,  equilibrium  requires  that  the  radial  and  circum- 
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\ 


ELASTIC 


Pb  >  Pb 

/ 


\ 


1 


Figure  1.  Outer  elastic  and  three  plastic  zones 


ferential  stress  components  are  equal,  that  is,  crf  =  (7»=  -pb-  If  the  material  is  elastic  and  in  a 
plane  strain  condition,  the  axial  stress  is  <T^~—2vp^  throughout,  where  v  is  Poisson’s  ratio 
(0  <  F  <  2).  Hence,  we  have  the  relationship  tr*  =  o-,  <  cTz.  With  increased  loading,  yielding  will 
occur  if 


where 


/  =  o-»  —  A/a-j  +  cr„  =  0 


1  +sin  » 
1  —  sin  ^ 


(1) 

(2) 


and  <7u>0  the  unconfirmed  crush  strength  tp  being  the  friction  angle.  Substitution  of  the  stress 
components  in  yield  condition  (1)  gives 


P.  =  Ph 


\-2Nv 


(3) 


For  Pb>0,  the  material  will  not  yield  if  2Nv>  1.  However,  if  2Nv<  1  the  material  will  yield 
throughout  when  equation  (3)  is  satisfied.  We  shall  therefore  consider  design  values  of  internal 
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pressure  p.  obeying  the  inequality 


P.<7  when  lNv<\  (4) 

1  —2Nv 

and  for  p„  so  far,  unrestricted  when  2Np  >  1. 

Once  the  design  internal  pressure  has  been  reached  it  is  held  constant  and  the  far-held 
pressure  is  increased.  While  the  material  remains  elastic  the  stress  and  strain  helds  about  a 
hole  of  radius  r  =  a  are 


a  a 

o'r=-Pb  +  (Pb-p.)-T-,  <F#=-Pb-(Pb-p.)-p-,  o-,  =  2i/pb  (5) 

a"  a" 

2Ge,=  -(l-2p)pb  +  (pb-p.)-p-.  2Cc«=-(l-2i')Pb-(pb-p.)^.  =  0  (6) 

where,  for  our  loading  history,  p,<Pb.  In  (5),  G  is  the  shear  modulus.  We  see  immediately 
from  (5)  that  at<cr,  and  at<cri  throughout.  Also,  <rr(oo)<crz,  but  at  the  hole  the  inequalities 
are  conditional,  that  is,  <T,(a)><Tt  if  p,<2vpb  and  tri<(T,(a)  if  p,<2j/pb.  It  can  readily  be 
shown  that  with  increasing  loading,  yielding  occurs  hrst  at  the  hole.  We  now  examine  the  effect 
these  inequalities  have  on  the  initial  yielding  at  the  hole. 

Initial  yielding 

Let  yielding  occur  at  /■  =  a  when  Pb  =  Pb.  If  p.  <  2ppb  when  yielding  occurs,  then  <  cr,  <  <r, 
at  r  =  a  and  the  yield  condition  is 


/  =  0-*  -  No-,  +  <7u  =  0  (7) 

Substitution  of  <T,(a)=  —  2pb  +  p,  and  cr,{a)-  — p»,  from  (5),  in  equation  (7)  leads  to  the  yield 
pressure 

Pb  =  J(N  +  l)p.-hcr„]  (8) 

Satisfaction  of  the  inequality  p»<2ppt,  now  requires 


W7u 

P.< - - - 

^  l-{N+l)p 


(N  +  1><1 


(9) 


but  places  no  restriction  on  positive  p.  when  (N  +  1)f>  1. 

If  2t'Pb<p,  <Pb  whc'  yielding  occurs,  then  cr«<(rr<crt  at  r  =  a  and  the  yield  condition  is 
equation  (1).  Substitution  of  cr,(a)=  -2pb  +  p,  and  az{a)=  —2ppb  in  equation  (1)  leads  to  the 
yield  pressure 


^*’”2(1 -Nr) 


Np<\ 


(10) 


Yield  condition  (1)  does  not  govern  yielding  if  Nr>l.  In  this  case,  pb  increases  until 
p,  <  2rpb  so  that  yielding  is  governed  by  condition  (7)  and  pb  is  given  by  formula  (8);  because 
(N+  I)r  >  1,  no  restriction  is  placed  on  p,.  For  Nr  <(N  +  l)r  <  1,  inequality  2rpb<p.  <Pb 
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and  the  yield  pressure  (10)  restrict  p.  to  the  ranges 


Wu 

— — - - <P.< 

,  0<Nv<\ 

(11) 

l-2Nv 

Wu 

l-(fv  +  l)»'*^^* 

(12) 

The  above  results  are  collected  in  Table  I. 


Table  I.  Initial  yielding  conditions 


We  present  here  the  analysis  of  Case  1  in  Table  I  for  Pb>Pb,  where  pb  is  given  by  (8).  In 
order  to  avoid  an  excessively  long  exposition,  Case  2  of  Table  I  will  be  presented  in  a 
forthcoming  paper. 

Initial  growth  of  plastic  zone 

As  the  far-field  pressure  is  increased  beyond  the  initial  yield  value  Pb,  an  annular  plastic 
zone  of  outer  radius  r  =  R  expands  from  the  hole  of  radius  r  =  a.  within  this  zone  yielding  is 
governed  by  condition  (7).  the  equation  of  equilibrium  throughout  the  material  is 

r^  +  cr,-o-,  =  0  (13) 

dr 

After  eliminating  from  (13)  by  means  of  yield  condition  (7)  it  is  found  that  the  solution  of 
the  resulting  equation  that  satisfies  a-r{a)=  -p.  is 

The  flow  rule  associated  with  yield  condition  (7)  is 

cj’’' =A— =-AN.  ei‘”  =  A— =  A.  ei‘”  =  A— =  0  (15) 

dcT,  d(Tt  dCTj 
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where  the  dot  may  be  interpreted  as  differentiation  with  respect  to  Pb-  From  flow  rule  (15)  it 
follows  that 

[eJ‘”+Ne‘,'”l  =  0.  e,‘'’*  +  N£S’'  =  g,(r).  ci'”  =  g2(r)  (16) 

which  means  that  the  strain  expressions  in  (16)  evaluated  at  a  fixed  radius  r  remain  constant  as 
loading  proceeds.  Initially,  when  the  elastic-plastic  radius  R  is  at  this  radius  r,  the  plastic 
strains  are  zero  so  that  gi(r)  =  0  and  g2(r)  =  0.  Hence, 


ci'’'  +  N£i‘’’  =  0.  £?’  =  0 


(17) 


Strain  increments  are  taken  as  the  sum  of  the  elastic  and  plastic  strain  increments,  that  is, 

£,  =  £r+£‘‘’',  =  +  =  +  (18) 

Now  c,-(e{*’  +  ci‘’*)  =  g(r).  and  when  R  =  r  we  have  E,=  ei*’  and  £^”'  =  0  so  that  g(r)=0. 
Using  the  same  argument  for  the  other  two  strain  components  leads  to 

£,  =  £:-'+£,'■”,  £.  =  £r  +  cr.  E,  =  £i'>  +  £l‘’’  (19) 

Plane  strain  requires  £i  =  0  and  the  flow  rule  leads  to  ci'’’  =  0,  as  shown  in  (17).  Thus  (19)  for 

the  sum  of  the  strains  completes  the  deduction  that 

£.  =  0,  £?’  =  0.  £i'’'  =  0  (20) 

Hooke’s  law  relating  the  elastic  strains  to  the  stresses  is 

Ee  =  <rt-  •'(<^0  +  iZi)  (2 1 ) 

Ee  i*'  =  -  »'(<T,  +  a,)  (22) 

fci*’ =<r,-»'(<7,  +  cr»)  (23) 

where  E  is  Young’s  modulus.  The  result  ci*’  =  0  substituted  in  (23)  gives 

cr,  =  i>{crr-¥cra)  (24) 

Elimination  of  and  tr,  from  (21)  and  (22)  by  usin^  yield  condition  (7)  and  result  (24)  leads 
to 

2G£‘*’  =  (1-{N  +  1)»/]o-,+  w„  (25) 

2Ge!r'=[N-(Ar  +  l)p]o-,-(l-p)<ru  (26) 

where  G  is  the  shear  modulus  and  cTr  is  given  by  (14). 

The  strains  are  related  to  the  radial  displacement  u  by 


du 

"'  =  d7 


U 

r 


(27) 


and  elimination  of  u  leads  to  the  compatibility  equation 

r^^+£»-£,  =  0  (28) 

dr 

By  using  in  turn  (19)  for  decomposing  the  strains,  (17)  to  eliminate  cj'’’,  (25)  and  (26)  to 
eliminate  the  elastic  strains  and  (14)  to  eliminate  <r,,  equation  (28)  becomes 


dr 


+  (N  +  1  )£  ?’  =  ^  •  (N  +  1  )((N  -  1  )p.  + 


(29) 
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The  solution  of  equation  (29)  satisfying  e«’'(/?)  =  0  is 


_<p>  —  _ 
£«  — 


l-i>  N+l 
2G  '  N 


and  e,  is  given  by  (17).  The  remaining  unknown  is  the  radius  R.  Using  R  and  p  instead  of  a 
and  p,  in  (5)  gives 

a-,(R)=-p,  a«(/?)=-2pb  +  p,  o-,(/?)= -2^pb  (31) 

If  yield  condition  (7)  still  holds  at  /•  =  /?  it  follows  that 

Pb  =  2[N  +  l)p  +  trul  (32) 

(33) 

/V  +  1 

Substituting  r  =  /?  in  formula  (14)  and  equating  to  formula  (33)  provides  continuity  of  the 
radial  stress  and  the  result 

2  (N-l)pb  +  tr„ 


2 

Va7  Af+1 


(N  -  1  )P,  +  CTu 


The  above  description  of  yielding  applies  with  continued  loading  as  long  as  <Ti<a,  in  the 
plastic  zone.  From  (24)  and  (7), 


a-,  — o-,  =  (l— (A^  +  l)i/]er,+  I'cr^ 


so  that 


(1)  when  (N  + !)»'>  1 

(2)  o-t<<T,  when  (N  +  1 )»' <  1  ifcr,<- 


l-(N  +  l)j/ 


In  the  first  case,  when  (A/  +  l)j/  >  1,  the  solution  is  complete;  this  case  has  been  treated  in  the 
literature.^ ^  In  the  second  case,  when  (A^+l)t'<l,  we  first  observe  that  in  the  plastic  zone 
min  O’,  =  <rr(R )  so  that  as  loading  proceeds  the  elastic-plastic  radius  R  attains  a  value  R  where 


o-,(«)  =  o-,(i?)=- 


l-(N  +  l)i/ 


o-,(«)  = 


(1  -t/)<T„ 

l-fN+l)*- 


The  magnitude  of  the  loading  when  R  =  R  is  Ph  given  by  (31)  and  (35),  that  is. 


2(1 -(N  + 1)^1 

and  the  radius  of  the  elastic-plastic  interface,  from  (34)  and  (37),  is  determined  by 

(1-2./)<7„ _ 

Va/  (I-(A/+1)i/J[(A/-1)p.  +  (7„I 
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Relationship  (38)  for  A  >s  independent  of  pt,  as  is  the  stress  held,  and  consequently  both 
remain  constant  during  further  loading. 

We  now  show  that  with  further  loading  three  plastic  zones  appear,  as  shown  in  Figure  1. 
Increased  loading — outer  plastic  zone 

It  is  postulated  that  with  increased  loading  (pb>Pb)  plastic  yielding  outside  the  fixed  radius 
A  and  adjacent  to  the  outer  infinite  elastic  zone  occupies  an  annular  region  A  <r<R  in  which 
(rt<a-,<a’z,  as  shown  in  Figure  1.  The  yield  condition  in  this  zone  is  therefore  equation  (1) 


and  the  associated  flow  rule  is 

£‘'■'  =  0.  ei‘’'  =  A.  (39) 

which,  for  the  same  reasoning  that  led  from  flow  rule  (15)  to  the  plastic  strain  relations  (17), 
provides  the  strain  relations 

£‘'’'+Ne‘/‘  =  0,  £j‘‘’  =  0  (40) 

Relations  (19),  the  plane  strain  condition  £<  =  0  and  relations  (40)  lead  to  the  strains 

£r=£r*\  - £1*'.  Cl  ~  0  (41) 

iV 

Hooke’s  law  (21)-(23),  expressions  (41)  and  yield  condition  (1)  give 

+  1  -  (AT  +  lfu]crr  =  (/V^  +  1  - 2Nu)2Ger  +  N(I^  +l)t>2G£»  +  (N-l )v(t^  (42) 

[N^+l-(N+  =  N{N  +  \)i>2Gb.  +  nHG£s  -  [1  -  (A/  +  Ot-ja-u  (43) 


[A/^  +  1  -  (N  +  l)\]cr»  =  N(N  +  l)u2Ge.  +  N^2G£e  -  [1  -  (A/  +  l)i-]c7-„  (43) 

Substitution  of  the  stresses  (42)  and  (43)  in  the  equation  of  equilibrium  (13)  and  use  of  the 


strain-displacement  relations  (27)  lead  to  the  displacement  equation 

,d*u  du  2  (l—2u)a-u 

^  dr^'^'  dr  °  “  N  2G 

where 

(44) 

a'  =  A/V(A/^+l-2A/i/) 

The  displacement  solutions  of  (44)  and  the  corresponding  strains  are 

(45) 

2Gu  =  A  ir”  +Air''’  +  Ar 

(46) 

2Gb,  =  ir"”'  -  aA^r~'‘~'  +  A 

(47) 

2C£,  =  /i,r"~'  +  A2r-"''  +  /\ 

where 

,  (l-2»')o-u 

(48) 

l-2Nu 

(49) 

If  the  elastic  zone  stresses  evaluated  at  r  =  R  are  substituted  in  yield  condition  (1)  we  find 
that 


o’r(R)= -2(1 -Nr)pb  +  cr„  (50) 

Hence,  the  circumferential  strain  at  r  =  R  is  given  by 

2G£,(R)=-2(N-l)..pb-(T„  (51) 
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The  radial  stress  and  circumferential  strain  at  r  =  from  the  plastic  zone  are  determined  by 
(42),  (47)  and  (48).  Thus,  to  provide  continuity  of  radial  stress  and  displacement  (or  f«)  at 
r  =  R  'nc  require,  according  to  (50)  and  (51), 

/4  t/g"-*  +  - n  - 2Ny)p„]  (52) 

\  —  2Nu 

[^+ (N  +  l)Fj  A  ,/g  (N  +  1  )I/] j/g 

2(l-Nu)[N^+\-{N+\)\],  . .  ,,,, 

- N(l^C) - [---(1-2N.)P.]  (53) 

It  is  reasonable  to  postulate  that  at  the  inner  radius  R  of  this  plastic  zone  we  have 
<r,(R)  =  ai(R)  so  that  the  yield  condition  may  be  written  as 

Na,(R)-a,(R)  =  o-.  (54) 

where  the  stresses  are  given  by  (42)  and  (43).  Substitution  of  (42)  and  (43)  in  (54)  and  use  of 
strain  solutions  (47)  and  (48)  with  r  =  R  lead  to  the  equation 

ia~  0[o ~(a  l][“(N  +  1)v]a2^'"“’ 

[N^+l-(N  +  l)^](l-2v)o-, 

a(l-2Ni')  ^  ’ 

After  solving  (52)  and  (53)  for  Ai  and  Az,  and  substituting  in  (55)  we  arrive  at  the  equation 
ia~  +  l)v]  [  j^-  (AT  +  l)t/](N -i)Nu  +  {N^+l-iN+l  )^f}(1  -  N./)] 

.  [ I (yv  +  1  )i/l(A^  -  1  )/V,/  -  +  1  - (A^  +  1  )V}(1  -  A^*^)]  (^) 

=  (N"+1-2Nv)(N'+1-(N  +  1)^]— (56) 

Equation  (56)  determines  the  radius  ratio  R/R.  Determination  of  R  requires  a  description  of 
the  state  in  the  plastic  zone  R  <r<R. 

Increased  loading — middle  plastic  zone 

It  is  postulated  that  with  increased  loading  {pb> Pb)cr,  =  cr^  throughout  the  plastic  zone 
R  <r<R  shown  in  Figure  1,  so  that  yielding  is  governed  by 

fi  =  crt—N<Ti  +  a-u  =  0  (57) 

/:  =  O'#  —  Nor,  +  cTu  =  0  (58) 

Equation  (57)  allows  the  elimination  of  o-#  from  the  equilibrium  equation  (13)  to  leave  a 
differential  equation  with  the  solution 


,  ,  N-I  _  ,  ,  N-I , 


71 


COMPRESSION  OF  A  MOHR-COULOMB  MEDIUM  AROUND  A  HOLE 


375 


which,  after  substituting  R  and  from  (38)  and  (35),  becomes 

Because  formulae  (60)  and  (14)  are  the  same,  the  radial  stress  tr,  is  given  by  the  same  formula 
throughout  the  combined  inner  and  middle  zones,  a<r<R.  Formula  (60)  gives  the  value 
cr,(R),  so  providing  continuity  of  this  stress  at  the  interface  r  =  R  leads  to  the  equation 

(N  +  Dt/j  [  j  (N  +  1  )«/}(N  -  1  )N«/  +  {Af ^  +  1  -  (yv  +  1  )^v}(l  -  yv./) j 

+  [^- (iV  + 1  )./j  1^  j  (N  + 1  )^  j (iV  - 1  )yv»/ -  + 1  -  (A/ + 1 ) Vk  1  - yvt.)]  ( 

(yv^+i-2yv./)[N^  +  i-(N  +  i)^i/]r„  ^ 

(61) 

The  radius  ratio  R/R  is  determined  by  (56),  so  equation  (61)  determines  the  radius  R  and 
hence  R  is  also  determined.  Now  that  R  is  known,  the  constants  A\  and  A2  (in  terms  of  the 
loading  parameter  pb)  are  the  determinable  solutions  of  (52)  and  (53).  Explicitly, 


A,  =  ;?'~“[1^-(N +  l)p](yV-l)Ni/  +  {N^+l -(A/+ !)>}(! -yvj/)] 
Ai  =  /?'*“[  j^+ (yv  +  1)./[(N  -  1)A/f  -  + 1  -  (N  +  l)^t'}(l  -  A/p)] 


a[q-u-(l  -2A/p)pb] 
A/"(l-2A/p) 

(^ 

a[cr.i-(l  -2iVp)pbl 
A/^l-2A/p) 


and  by  (46),  (47)  and  (48)  the  displacements  and  strains  in  the  outer  plastic  zone  are  known. 
Finally,  by  (42),  (43)  and  (1)  the  stresses  in  the  outer  plastic  zone  are  known. 

The  flow  rule  associated  with  the  yield  condition  equations  (57)  and  (58)  is 


,-(p)  _  1 
£  r  ”  1 

—  +  A  2 — “  ”AfA  1 
acr,  da. 

(64) 

^fP)  _  1 

df\  a/2 

— +  A2— =  A,+A2 

da»  daa 

(65) 

,;tP)  _  1 
£i  —  A  1 

da,  da. 

(66) 

[£j'”+A/ey”  +  £i‘” 

]=0,  £:'’'+A/£'/’+£'‘’'=g(0 

(67) 

so  that 


Let  the  plastic  strains  he  fj”’,  e»’’,  el”*  at  r  when  the  load  Pb  has  placed  the  plastic-plastic 
inte'rface  of  the  middle  and  outer  zones  at  r  Then,  when  R  =  r, 

g(r)=£:'”  +  A/fr  +  4'’‘  (68) 

But  in  the  outer  plastic  zone  the  plastic  strain  relations  (40)  hold  so  (68)  reduces  to  g(r)  =  0. 
Hence  (67)  becomes 

£['”+Ncir’+ci'’'=0  (69) 
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The  plane  strain  condition  is  e,  =  ti'’  +  Cz'’*  =  0  which,  when  combined  with  (69)  gives 

=  +  +  (70) 

Stress  relationships  (57)  and  cr,  =  o-i  reduces  Hooke's  law  (21),  (22)  and  (23),  to 

Ee  '“  =  [  1  -  (N  +  1  )v)o-z  +  t'o-u  (71) 

£e‘/’  =  (N-2./)o-,-o-„  (72) 


where  <r,  is  given  by  (60).  Substitution  of  the  elastic  strains  (71)  and  (72)  in  the  radial  strain 
expression  (70)  followed  by  substitution  of  the  resulting  e,  in  the  compatability  equation  (28) 
leads  to  the  equation 


dr 


N^  +  2-2(2N  +  \)u, 


(N-l)(l  +  v) 


2G—  (r  £,)= - — — ,,,,  , —  [(N  -  1  )pa  +  (Tu)  -^frrr  ---  .  ■  r  (73) 


(N  +  2)(l-2v)o-„ 


(N-l)(l  +  v) 


Integrating  (73)  from  r  to  R  gives 


.  Ar^  +  2-2(2Af  +  l)r 
2Ge.  =  2Ge,(R).  (-)  + 


(Af  +  2)(l-2v)o-u 
(N-lXN  +  l)(l  +  v)L\r  J 


(74) 


where  2Ge»(£)  is  determined  by  (48).  From  (70),  (71)  and  (72)  the  radial  strain  is  determined 
by 


,,  ,  A/^  +  2-2(2iV+l)v  W-2v 

2c7cr  =  -2GeJ\  + - T" - (T, — : - a-^ 


1  +  ) 


1  +  V 


(75) 


Increased  loading — inner  plastic  zone 

During  the  increase  in  loading  (pb>Pb)  the  stress  field  in  the  inner  plastic  zone,  a  <r<R 
remains  constant.  Formula  (14)  gives  <t,.  Formulae  for  <7»  and  o-j  then  follow  from  the  yield 
condition  (7)  and  the  result  (24).  The  differential  equation  (29)  for  the  plastic  strain  eg"’  still 
holds  with  increased  loading  but  now  we  integrate  from  rtoR  instead  to  r  to  R.  Instead  of  (30) 
we  obtain 


2Ge 


(p)  _ 


=  2G£‘,'”(i? 


(N+l)(l-v) 

2N 


,  Dn  N-lr  ,  i5,  Nvl  ,N-1, 

i(A/-np.+<7.i(^  [(7)  -ij)  ] 


(76) 


The  other  piastic  strain  components  are  still  given  by  (17).  The  elastic  strain  components 
remain  the  same  because  the  stress  field  is  unchanged,  so  the  elastic  strains  are  given  by  (20), 
(25)  and  (26).  In  (76),  e^f^{R)  is  found  by  setting  r  =  £  in  (74)  to  give  eg(f? ),  r  =  R  in  (72)  to 
give  £t\R)  and  applying  strain  summation  (19);  note  that  for  (72)  (t,{R)  is  given  by  (35). 


Numerical  example  (Case  /,  Table  1) 

A  numerical  example  is  presented  to  illustrate  a  comparison  between  the  stress  and  strain 
fields  obtained  by  the  analysis  and  the  NONSAP  finite  element  code  augmented  by  the 
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Mohr-Coulomb  yield  condition.  The  data  are 

1^=  1/5 

N=  3(^3  =  30°) 
tr„=  2000  Ib/in^ 
pa  =  500  Ib/in* 


Poisson’s  ratio 
Friction  angle  function 
Unconfined  crush  strength 
Internal  pressure 


We  can  load  initially  with  p,  ~  pb  until  p,  =  500  Ib/in^  without  yielding  because  2Nf  =  6/5  > 
1;  see  condition  (4).  The  chosen  numerical  values  result  in  (N  +  l)t3  =  4/5  <  1  and 
(Af  +  l)j/)  =  2000  Ib/in^,  so  that  inequalities  (9)  are  satisfied.  Hence  pi<2vpt,  and  the  initial 
yield  pressure  is  determined  by  formula  (8),  which  gives  Pb  =  2000  Ib/in'  {2vpb  =  800  Ib/in^). 
The  value  of  the  far-field  loading  when  R  =  R,  the  radius  of  the  elastic-piastic  interface  when 
a-,{R)-cri(R),  is  Pb  =  5000  Ib/in^  according  to  formula  (37);  the  stresses  at  r=  R  are  cr,{R)  = 
ai(,R)=  -2000  Ib/in^  and  o-*(R)  =  —8000  Ib/in^,  accoiding  to  formulae  (35)  and  (36).  Formula 
(34)gives  R  =  'J(2)a.  Another  result  that  can  readily  be  calculated  is  the  circumferential  strain  at 
the  hole  when  pb  =  Pb-  According  to  formulae  (26)  and  (30),  2C£’5’(fl)  = -2700  Ib/in'  and 
2Gf /’(a)  =  - 1 1 ,200  Ib/in^  so  that  2G£g{a)=  13,900  Ib/in'.  If,  for  simplicity,  G  =  10^  Ib/in', 
£9(0)  =  0-70  per  cent  when  Pb  =  Pb  =  5000  Ib/in',  that  is,  the  tunnel  diameter  has  been  reduced 
by  0-70  per  cent. 
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RADIUS  RATIO,  r/a 


(94  NODES  AND  23  ELEMENTS) 

Figure  4.  NONSAP  axisymmelric  finite  element  model 
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For  loading  above  transcendental  equations  (56)  and  (6 1 )  have  to  be  solved  to  find  R  and 
R  before  stress  and  strain  formulae  can  be  evaluated.  However,  in  many  cases  the  constant  a 
of  formula  (45)  is  approximately  unity  and  rapid  estimations  can  be  made. 

Figures  2  and  3  afford  a  comparison  of  the  stress  and  strain  fields  for  a  far-field  pressure  of 
Pb  =  6000  Ib/in';  for  the  strain  calculations  the  modulus  of  rigidity  was  taken  as  C  = 
lO^Ib/in^.  Figure  4  shows  the  finite  element  model  used  in  the  NONSAP  calculations.  It  is 
clear  from  the  excellent  agreement  that  for  problems  of  this  type  the  Mohr-Coulomb  yield 
condition  was  correctly  integrated  into  the  NONSAP  finite  element  code. 
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FLORENCE  AND  SCHWER  CASE  II  SOLUTION 

The  following  paper  (Reference  4)  was  prepared  under  a  DNA  contract.  It 
is  reproduced  here  in  its  entirety  with  the  concurrence  of  the  senior  author. 
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SUMMARY 


An  analytical  solution  is  presented  for  the  stress  and  strain  fields 
in  a  Mohr-Coulomb  material  in  plane  strain  around  a  circular  hole  when  it 
is  compressed  by  an  axisymmetric  far-field  pressure.  Several  solutions 
arise  involving  one  to  three  plastic  zones  depending  on  the  material 
properties  and  the  loading.  A  solution  Involving  three  plastic  zones  was 
presented  in  reference  1.  The  solution  presented  here  involves  two 
plastic  zones.  The  solution  was  used  to  further  validate  the  functioning 
of  the  Mohr-Coulomb  yield  condition  that  was  added  to  the  NONSAP  finite 
element  code.  Stress  and  strain  field  comparisons  are  made. 
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INTRODUCTION 


The  purpose  of  this  paper  is  to  provide  an  analytical  solution  of  a  rock 
mechanics  problem  that  can  be  used  to  validate  the  predictions  of  finite  ele¬ 
ment  codes.  The  problem,  which  is  of  interest  in  itself,  is  the  determination 
of  the  stress  and  strain  fields  surrounding  a  tunnel  in  a  Mohr-Coulomb  rock 
subjected  far  from  the  tunnel  to  an  axisynroetric  pressure.  Plane  strain  con¬ 
ditions  prevail  and  an  internal  pressure  is  allowed  to  develop  at  the  hole 
boundary.  Such  a  pressure  may  be  considered  as  arising  from  a  metal  liner 
at  yield  or  from  a  crushable  back-packing  material  between  the  liner  and  the 
rock.  To  simplify  the  analysis  and  to  provide  a  defined  loading  history 
this  internal  pressure  is  assumed  to  be  equal  to  the  far-field  pressure 
as  the  far-field  pressure  is  increased  (loading)  until  the  Internal 
pressure  reaches  its  final  value,  whereupon  it  remains  constant  with, 
subsequent  loading. 

The  rock  behavior  is  elastic-perfectly  plastic  obeying  the  Mohr-Coulomb 
yield  condition  and  associated  flow  rule;  hence,  dilatation  is  included.  It 
will  be  seen  that  several  cases  can  arise  that  depend  on  the  relative  magni¬ 
tudes  of  the  principal  stresses  and  these  relative  magnitudes  depend  on  the 
values  of  Poisson's  ratio  and  the  friction  angle.  One  of  these  cases  is 
treated  in  reference  1.  It  involves  two  adjacent  faces  and  the  common 
edge  of  the  yield  pyramid,  and  three  corresponding  plastic  zones.  The  case 
treated  here  involves  one  face  and  an  edge,  and  two  corresponding  plastic 
zones,  as  shown  in  Figure  1.  In  this  work  the  finite-element  code  employed 
is  NONSAP^  to  which  the  Mohr-Coulomb  yield  condition  presented  in  reference  3 
has  been  added.  Comparisons  are  made  of  stress  and  strain  fields  obtained  by 
analysis  and  code. 
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The  motivation  for  the  code  validation  analysis  was  to  provide 
support  for  the  experimental  and  theoretical  investigations  of  Lindberg 
and  Kennedy**  and  Lindberg  and  Senseny®  into  the  response  of  buried 
cylindrical  structures  to  quasi-static  loading, 

INITIAL  LOADING.  OUTER  ELASTIC  ZONE 


Let  p  be  the  pressure  acting  on  the  hole  boundary  of  radius  a  and  let 
a 

p  be  the  pressure  at  infinity  (far-field).  If  p  ”  Pu  on  initial  loading, 

w  S  D 

equilibrium  requires  that  the  radial  and  circumferential  stress  components 

are  equal,  that  is,  =  Cg  =  -  p^.  If  the  material  is  elastic  and  in  a 

plane  strain  condition,  the  axial  stress  is  o^  =  throughout,  where 

V  is  Poisson's  ratio  (0  <  v  <  1/2).  Hence,  we  have  the  relationship 

O’-,  *  CT  <  a  .  With  increased  loading  yielding  will  occur  if 
o  r  2 


f*a.  -N<?  +  a  =  0 

0  2  u 

(1) 

where 

1  +  since 

N  =*  , - : — * 

1  -  sincp 

(2) 

and  a  >  0  is  the  unconfined  crush  strength,  ‘P  being  the  friction  angle.  Substi¬ 

tution  of  the  stress  components  in  yield  condition  (1)  gives 


P 


a 


(3) 


For  p,  >  0,  the  material  will  not  yield  if  2Nv  >  1.  However,  if  2Nv  <  1  the 

D 

material  will  yield  throughout  when  Eq.  (3)  is  satisfied.  We  shall  there¬ 
fore  consider  design  values  of  internal  pressure  p^  obeying  the  inequality 

a 

P  when  2Nv  <  1  (4) 

a  1  -  2Nv 

and  for  p  ,  so  far,  unrestricted  when  2Nv  >  1. 
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Once  the  design  internal  pressure  has  been  reached  it  is  held  constant 
and  the  far-field  pressure  is  increased.  While  the  material  remains  elastic 
the  stress  and  strain  fields  about  a  hole  of  radius  r  *  a  are 


(7^  =  +  (Pb  ■ 


■  Pb  ‘  ^Pb  -  Pa^~ 


'r  “  -  (1  '  ^''^Pb  ^Pb  •  Pa^^ 


(1  -  2.v)pj^  -  (p^^  -  p^)—  =  0  (6) 


where,  for  our  loading  history,  p^  <  p^.  In  (5),  G  is  the  shear  modulus.  We 

see  immediately  from  (5)  that  <  a  and  a.  <  o  throughout.  Also, 

or  0  2 

a  (»)  <  a  ,  but  at  the  hole  the  inequalities  are  conditional,  that  is, 
r  z 

Q  (a)  <  O  if  p  >  21^  and  O  <  a  (a)  if  p  <  2Up  .  It  can  readily  be 
rzabzr  a  b  ^ 

shown  that  with  increasing  loading  yielding  occurs  first  at  the  hole. 

We  now  examine  the  effect  these  inequalities  have  on  the  initial  yielding 


at  the  hole. 


INITIAL  YIELDING 


Let  yielding  occur  at  r  =  a  when  p^  =  p^.  If  p^  <  2vp^  when  yielding 
occurs,  then  Og  <  o^  <  at  r  =  a  and  the  yield  condition  is 


f  =  CT-.-Na  +  CT  =  0 
9  r  u 


Substitution  of  a-(a)  =  -  2  p.  +  p^  and  a  (a)  =  -  p  ,  from  (5),  in  Eq.  (7) 

u  D  o  r  o 


leads  to  the  yield  pressure 


Pb  “  2^^^  +  ^>Pa  ^  '"J 


Satisfaction  of  the  inequality  p^  <  2vpj^  now  requires 


1  -  (N  +  l)v 


(N  +  1)  V  <  1 


but  places  no  restriction  on  positive  p^  when  (N  +  l)v  >  1 
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If  2vpj^  ^  yielding  occurs,  then  Og  <  <  <7^  at  r  *  a  and 

the  yield  condition  is  Eq.  (1)  .  Substitution  of  CToCa)  =  -2p  +  p  and 

w  b  a 

(7  (a)  =  -21<'p  in  Eq.  (1)  leads  to  the  yield  pressure 
z  b 


p  +  o 
a  u 

2(1  -  Nv) 


Nv  <  1 


(10) 


Yield  condition  (1)  does  not  govern  yielding  if  Nv  >  1.  In  this  case, p^  in¬ 
creases  until  p  <  2vp.  so  chat  yielding  is  governed  by  condition  (7)  and  p, 
is  given  by  formula  (8);  because  (N  +  l)v  >  1,  no  restriction  is  placed  on  p 

3 

For  Nv  <  (N  +  l)v  <  1,  inequality  Ivp^^  ^a  ^b  yield  pressure  (10) 

restrict  p^  to  the  ranges 


va 


1  -  (N  +  l)v 


1  -  2Nv 


0  <  Nv 


<i 


(11) 


va 


1  -  (N  +  l)v 


^  <  Nv  < 


(N  +  l)v  <  1 


(12) 


The  analysis  of  case  1  in  Table  1  is  contained  in  reference  1.  Here, 
we  present  the  analysis  of  case  2  for  <  p^,  where  is  given  by  (10). 

INITIAL  GROWTH  OF  PLASTIC  ZONE 

As  the  far-field  pressure  is  increased  beyond  Che  initial  yield  value 

p  an  annular  plastic  zone  of  outer  radius  r  =  R  expands  from  the  hole  of 
b 

radius  r  =  a.  Within  this  zone  yielding  is  governed  by  condition  (1) 
and  the  associated  flow  rule 


(13) 


where  Che  docs  may  be  interpreted  as  differentiation  with  respect  to  Che 

loading  p,  .  From  flow  rule  (13)  it  follows  that 
b 
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INITIAL  YIELDING  CONDITIONS 
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(14) 


=  g  (r) 

y  z  1 


-  g  (r) 
t  2 


which  means  that  the  strain  expressions  in  (14)  evaluated  at  a  fixed  radius  r 

remain  constant  as  loading  proceeds.  Initially,  when  the  elastic-plastic 

radius  R  is  at  this  radius  r,  the  plastic  strains  are  zero,  so  that 

g  (r)  =  0  and  g  (r)  =  0.  Hence  the  strain  expressions  in  (14)  become 
1  z 


I'o'  +  •  0 


Ne'"'  +£'"'=  0  e 
6  z  r 


(P) 


(15) 


Strain  Increments  are  taken  as  the  sum  of  the  elastic  and  plastic 
strain  increments,  that  is. 


r  r  ^  r  ®e  ®e 


e  =  +  c^P^ 

z  z  z 


(16) 


Now  e  -  (e^®^  +  e^P^)  *  2(r).  and  when  R  *  r  we  have  e  =  and 

r  r  P  '  S'*-  ^  I  r  r 

£^^^  "0,  so  that  g(r)  *  0.  Using  the  same  argument  for  the  other  two 


strain  components  leads  to 


®r  ®r  ®e  ®0 


e  *  c 
z  z 


(e)  +  g(p) 


(17) 


Plane  strain  requires  £  =  0,  so  (15)  and  (17)  give  the  plastic  and  total 

z 

strains  as 


(a-)  ‘'P^ 

E^P^  =0  N£p,  =  £ 

r  9  z 


e^P>  = 
z  z 


and 


(e) 


,(e)  ^  1  /e) 


N  z 


(18) 
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Hooke's  law  relating  the  elastic  strains  to  the  stresses  is 


EC-'  “ 


E  =  CT  -  ■*■  o'e^ 


where  E  is  Young's  modulus.  Elimination  of  from  (19),  (20),  and  (21) 
by  using  the  yield  condition  (1) ,  and  substitution  of  the  resulting  elastic 
strains  in  (18)  gives  the  stress-strain  relationships 

[N^  +  1  -  (N  +  D^vla^  =  (N^  +  1  -  2Nv)2Ge^  +  N(N  +  l)v  2Geg  +  (N  -  D 

[N^  +  1  -  (N  +  U^vlOg  =  N(N  +  l)v  2Ge^  +  N^2Geg  -  [1  -  (N  +  l)v]a^ 

Substitution  of  these  stresses  in  the  equation  of  equilibrium 

da 

r  ^  +  a  -  a  *  0 

dr  r  o 

and  substitution  of  the  strain-displacement  relations 


®r  dr 


u 

=  ■; 


leads  to  the  radial  displacement  equation 

2  .  ^  ^2  (1  -  2l0a 

2  d  u  ,  ^  ^  cy2  -  £ _ y 

dr 

where 


=•  N^/(N^  +  L  -  2Nv) 
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The  displacement  solution  of  (26)  and  the  corresponding  strains  are 


a  -a 

2Gu  =  A^r  +  A^r  +  Ar 


or-l  ^  -Of-l 

2G€  *  oA.r  -  aA_r  +  A 

r  1  1 


IGCg  =  A^r“-^  +  A^r'^'^  +  A 


where 


(1  -  2v)ct 
_ u 

1  -  2Nv 


In  the  elastic  zones  the  stresses  and  strains  are  given  by  (5)  and 

(6)  with  a  =  R  and  p  =  -0  (R) ,  Thus,  yield  condition  (1)  applied  at  the 

d  IT 

elastic-plastic  interface  r  =  R  gives 


a  (R)  =■  -  2(1  -  Nv)p^  +  o 

r  u 


2G€g(R)  *  -  2(N  -  l)vp^  - 


Continuity  of  displacement  at  r  =  R  is  assured  by  equating  (33)  to  (30) 
and  continuity  of  radial  stress  at  r  =  R  is  assured  by  equating  (32)  to 
the  stress  in  (22)  after  eliminating  the  strains  with  (29)  and  (30).  These 
operations  provide  two  algebraic  equations  for  A^  and  A^  with  the  solution 


t)^(l  -  2SV) 


BoR  - 

N^(l  -  2Nv) 
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where 


(36) 

(37) 

The  radius  R  of  the  elastic-plastic  interface  is  found  by  satisfying 
the  pressure  boundary  condition  at  the  hole.  Hence,  (22)  with  a  =  -p  , 
and  with  the  strains  (29)  and  (30)  at  r  =  a  having  and  given  by 
(34)  and  (35),  provides  the  equation 


®1  "  [a  ■  +  [n^  +  1  -  (N  +  D^v]  (1  -  Nv) 

®2  “  [f  -  [n^  +  1  -  (N  +  D^v]  (1  -  Nv) 


(N  +  l)v 


1-a 


-f[» 


^  +  1  -  (N  +  D^v 


-  (1  -  2Nv)p^ 

-  (1  -  2Nv)p^^ 


(38) 


for  the  radius  R,  and  being  given  by  (36)  and  (37) . 

The  above  description  of  yielding  applies  with  continued  loading  as 
long  as  in  the  plastic  zone.  A  numerical  study  of  the  stress 

difference  0  -  a  was  performed  for  values  of  N  and  v  that  satisfied  the 

Z  IT 

inequality  Nv  <  (N  +  l)v  <  1  and  it  was  found  that  o  -a  decreases  as  p 

z  r  b 

increases  and  becomes  zero  first  at  the  hole,  r  =  a.  Explicit  expressions 
describing  the  stress  field  are 


Co  =  N  f  - 
r  La 


+  (N  +  l)v 


^  -  B  [a .  („  +  i)v] 


00 


(39) 


l+a  1-2NV 
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“e  -  ON  [I  +  <»  +  »'■]  -fe  +  [j  -  +  »'■]  S 

^  r  r 

C0^.a[S+  (H+l)v]^  +  <.[f-  (»4.i,v]^-^ 


A.  C2vO 
2  _ u 

1+a  "  1-2NV 


where 


C.  -  +  1  -  (N  +  D^v 


If  we  set  o  (a)  «  a  (a)  =  -p  in  (39)  and  (Al)  we  obtain  two  equations 

L  2  d 

for  the  loading  and  the  corresponding  elastic-plastic  radius  R*  that 
may  be  put  in  the  form 


.  h  [f  *  ^  -  (l  - 

*  ^  B/[f  -  (N  +  l)vj  +  ijo^  .  +  lj(l  .  2»v)p/] 


\  -  <1  -  2Nv)p;] 


^  -  f  *  i) 


ABiB^a^  (  a  “  +  l)v 


+  (N  +  1)V  ^ 


INCREASED  LOADING 

It  is  postulated  that  with  increased  loading  (p.  >  v')  a  radius 

b  b 

R  >  a  exists  where  0^(R)  =•  0^(R).  In  the  outer  annular  plastic  zone 

R  <  r  <  R  we  still  have  Oq  <  O  <  a  (Figure  1)  and  the  stress  and  strain 

0  r  z 

fields  are  given  by  (39)  -  (A2)  and  (29)  -  (31)  but  the  radius  R  >  R*  is 

no  longer  given  by  (38).  In  fact,  satisfying  the  condition  o  (R)  “  o  (R) 

^  z 
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by  equating  (39)  and  (41)  with  r  =  tv  leads  to 


for  the  radius  ratio  R/R.  The  individual  values  depend  on  the  fields 
in  the  zone  a  <  r  <  R. 


It  is  also  postulated  that  with  increased  loading  0^  =  in  the 
plastic  zone  a  <  r  <  R,  as  shown  in  Figure  1,  so  that  yielding  is  governed 
by 

f^  =  Oq  -  N0^  +  CTu  =  0  (^6) 

f^  “  O-  -  Na  +0  =>0  (47) 

2  9  z  u 


Equation  (46)  allows  the  elimination  of  from  the  equilibrium 

y 

equation  (24)  to  leave  a  differential  equation  with  the  solution 


0 

r 


(48) 


which  satisfies  the  pressure  boundary  condition  at  the  hole.  Continuity 
of  the  radial  stress  at  the  interface  r  =  R  is  assured  by  setting  r  =  R 
in  (48)  and  (39)  and  equating  the  resulting  expressions  for  0^(R).  This 
operation  leads  to  the  following  equation  for  R  in  terms  of  the  ratio 
R/R,  which  is  the  solution  of  (45), 
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(N  -  Dp  +0 
a  u 


N(1  -  2v)a^ 
1  -  2NV 


<'?  +  +  1,„  1  Bilf -12  .  (»  H.  l)v{ 


1\R 


a 


a(N  -  l)[a  -  (1  -  2Nv)p^] 

U  D 


CN(1  -  2Nv) 


(49) 


The  flow  rule  associated  with  the  yield  condition  given  by  (46)  and 
(47)  is 


r  1  ^*2 


(50) 


/  X  ^^2 


(51) 


:(p)  .  X  ^  *  -  NX. 

^1  do  2  do 


(52) 


so  that 


r  9  z 


+  de'l’’  +  e'"'  =  8(r)  <53) 

r  u  z 


Let  the  plastic  strains  be  at  radius  r  when  the  load 

r  u  z 

has  placed  the  plastic-plastic  interface  at  r.  Then,  when  R  =  r 


e^P>  +  NelP^  +  i^P^ 


g(r)  -  e;*"  +  Ne^*'  + 


(54) 


But  in  the  outer  plastic  zone  the  plastic  strain  relationships  (15) 
hold  so  (54)  reduces  to  g(r)  =  0.  Hence  53)  becomes 


:^P^  +  Ne^P^  +  c^P^ 

r  y  2 


(55) 
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The  plane  strain  condition  is  e  *  »  0  which,  when 

z  2  z 

combined  with  (55)  gives 


S  -  -S  *  "s'*'  * 


The  yield  condition  (A6)  and  a  =  a  reduces  Hooke's  law  (19),  (20), 

2  ? 


and  (21)  to 


Ee^®^  =■  Ee^®^  =  ( 1  -  (N  +  l)v]a  +  va 
r  z  '  ^  T  u 


Eel®^  =  (N  -  2v)a  -  a  (; 

0  r  u 

where  is  given  by  (48) .  Substitution  of  the  elastic  strains  (57)  and 
(58)  in  the  radial  strain  expression  (56)  followed  by  substitution  of  the 


resulting  in  the  compatibility  equation 


leads  to  the  equation 


d7 


N~  4-  2  -  2(2N  l)v 
(N  -  1)(1  +  v) 


2N-1 

1  Vt" 

a  u  N-I 


(N  +  2)  (1  -  2v)ct 
_ u 

(N  -  1)(1  v) 


Integrating  (60)  from  r  to  R  gives 


.  xN+1 


2G€g  -  2G€g(R) 


.  +  2  -  2(2N  +  l)v . .  ,r/R\^” 

2N(N  -L)a  -rv) 


(N  +  2)(1  -  2v)CT 
(N  -  L)(N  +  1)(1 


-  2v)<7  r/p\N+’-  1 


9? 


(62) 


where  ZGe^CR)  is  determined  by  (30).  Explicitly, 

W 


2Geg(R) 


a  [a^  -  (1  -  2Nv)p^] 
N^(l  -  2Nv) 


(1  -  2v)a 
_ u 

1  -  2NV 


Substitution  in  (56)  of  the  elastic  strains  (57)  and  (58) ,  with 
given  by  (48) ,  results  in  the  radial  strain  formula 


2Ge^  -  -2GNeg  -  t(H  -  1)P,  (J) 


N-1  (N  +  2)(1  -  2v)a 


(N  -  1)(1  +  V) 


During  the  analysis  we  postulated  that  a  third  plastic  zone  could 

exist  next  to  the  hole  in  which  a„  <  O  <  o  ,  as  in  case  1  of  Table  1.^ 

9  2  r' 

The  result  of  the  analysis  showed  that  0  “  CJ  in  the  postulated  third 

z  t 

zone  so  that  only  two  plastic  zones  exist  in  case  2. 


NUMERICAL  EXAMPLE 


A  numerical  example  is  presented  to  illustrate  a  comparison  between 
the  stress  and  strain  fields  obtained  by  the  analysis  and  the  NONSAP 
finite  element  code  augmented  by  the  Mohr-Coulomb  yield  condition.  The 
data  are 


Poisson's  ratio 
Friction  angle  function 
Unconfined  crush  strength 
Internal  pressure 


V  -  1/5 

N  -  2  (<1^=  19.5°) 

0  =»  3000  psi 

u 

p  =  2000  psi 
3 


The  values  of  Poisson's  ratio  and  friction  angle  are  such  that 

2Nv  “  4/5  <  1  so  the  yielding  with  p  =  p,  does  not  occur  until 

a  D 
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p  *  P,  *  15,000  psi,  according  to  (3).  Hence  we  can  load  with  equal  p 

3d  3 

and  p  until  p  =  P,  =*  2000  psi.  Inequality  (11)  is  1500  psi  <  p  <  15,000  psi. 

D  3  D  3 

If  we  now  fix  p^  and  increase  p^  yielding  occurs  when  reaches  *  4167  psi, 

according  to  (10) .  As  loading  is  increased  a  single  plastic  zone  spreads 

until  ®  “  4529  psi  and  R  /a  =  1.10,  according  to  (43)  and  (44). 

Further  loading  creates  two  plastic  zones  and  when  p^  *  6000  psi,  we  obtain 

R/R  *  1.31  from  (45),  R/a  =  1.19  from  (45)  and  (49),  and  hence  R/a  =  1.55. 

The  tunnel  closure  at  this  loading  is  0.86%,  according  to  (61). 

In  many  cases  the  constant  a  of  formula  (27)  is  approximately  unity. 

This  value  allows  all  unknowns  to  be  determined  explicitly  so  that  rapid 

estimations  can  be  made. 

Figures  2  and  3  provide  a  comparison  of  the  stress  and  strain  fields 
as  obtained  by  the  analysis  and  by  the  NONSAP  finite-element  code.  Figure  4 
shows  Che  finite  element  model  used  in  the  NONSAP  calculations.  The  far- 
field  pressure  is  p^  =  6000  psi  and  a  modulus  of  rigidity  of  G  =  10®  psi  was 
taken  for  the  strain  calculations.  It  is  clear  from  the  excellent  agreement 
that  for  problems  of  this  type  the  Mohr-Coulomb  yield  condition  was  correctly 
integrated  into  the  NONSAP  finite-element  code.  The  same  conclusion  was 
reached  in  reference  1  for  case  1  of  Table  1. 
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ATTN:  M  JOHNSON 

INSTITUTE  FOR  DEFENSE  ANALYSES 
AHN:  CLASSIFIED  LIBRARY 
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DNA-TR-91-38  (DL  CONTINUED) 


KAMAN  SCIENCES  CORP 

ATTN:  LIBRARY  B  KINSLOW 
ATTN:  RICHARD  KEEFFE 

KAMAN  SCIENCES  CORP 
ATTN:  DASIAC 

KAMAN  SCIENCES  CORPORATION 
ATTN:  DASIAC 

LACHEL  AND  ASSOCIATES.  INC 
ATTN:  J  BECK 

NTS  ENGINEERING 

ATTN:  R  KENNEDY 
ATTN:  S  SHORT 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  D  WILSON 
ATTN:  HBRODE 

PACIFIC-SIERRA  RESEARCH  CORP 
ATTN:  DGORMLEY 

R  &  D  ASSOCIATES 

ATTN:  LIBRARY 

R  &  D  ASSOCIATES 

ATTN:  T  EDWARDS 

R  4  D  ASSOCIATES 

ATTN:  D  PIEPENBURG  C/0  TERRI  BOZARTH 

RAND  CORP 

ATTN:  B  BENNETT 

S-CUBED 

ATTN:  KDPYATT.JR 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  DR  M  MCKAY 


ATTN:  H  PRATT 

ATTN:  TECHNICAL  REPORT  SYSTEM 

SCIENCE  APPLICATIONS  INTL  CORP 
AHN:  WLAYSON 

SOUTHWEST  RESEARCH  INSTITUTE 
ATTN:  B THACKER 

SRI  INTERNATIONAL 
ATTN:  JGRAN 

THE  PBMK  TEAM 

ATTN:  DOCUMENT  CONTROL 

THE  TITAN  CORPORATION 

ATTN:  KKREYENHAGEN 
ATTN:  LIBRARY 
ATTN:  MITO 
ATTN:  S  SCHUSTER 

TRW  INC 

ATTN:  BBALACHANDRA 

TRW  SPACE  4  DEFENSE  SECTOR 
ATTN:  PNMATHUR 

TRW  SPACE  4  DEFENSE  SECTOR 

ATTN:  OUT6/W  WAMPLER 

UTD.  INC 

ATTN:  E  FOSTER 

WEIDLINGER  ASSOC,  INC 
ATTN;  FWONG 
ATTN:  H  LEVINE 

WEIDLINGER  ASSOCIATES,  INC 
ATTN:  I  SANDLER 
ATTN:  M  BARON 
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